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CHAPTER  I 


INTRODUCTION 


Although  one  cannot  expect  exact  solutions  for  a given  quantum 
mechanical  eigenvalue  problem,  it  is  often  possible  to  find  upper  and 
lower  bounds  to  the  eigenvalues.  Upper  bounds  can  be  obtained  satis- 
factorily by  the  well-known  variational  calculations,  for  example, 

[1  2l 

Rayleigh-Ritz  procedures  ’ ; however,  in  order  to  estimate  the 

accuracy,  the  knowledge  of  lower  bounds  is  essential. 

The  problem  of  finding  lower  bounds  has  been  approached  in 
many  different  ways.  Among  those,  the  intermediate  problems  origi- 

[31  [4] 

nated  by  A.  Weinstein  and  extended  and  developed  by  Aronszajn 

have  been  widely  used.  Many  fundamental  problems  in  quantum  theory  have 

[5l 

been  studied  by  Bazley  and  Fox  using  these  intermediate  problems. 

Recently  P.  0.  LBwdin  has  developed  a new  method^^-’'^^  for  the 

study  of  lower  bounds  to  energy  eigenvalues  by  the  use  of  the  parti- 
r Q 1 

tioning  technique^  . In  the  partitioning  technique,  the  first  order 
iteration  procedure  ^ ~ brackets  at  least  one  true  eigen- 

value E.  By  means  of  an  "inner  projection"^^^  on  a finite  manifold 
in  the  Hilbert  Space,  one  can  evaluate  lower  bounds  to  any  desired 
accuracy.  For  a Hamiltonian  of  the  form  + V,  where  V is 

positive  definite,  one  may  consider  the  intermediate  Hamiltonians: 


^ (n)  = ^ 


1 


where  Q is  a Hermitian  projection  operator  for  an  arbitrary  subspace 
of  order  n,  which  for  n-^oo  converges  towards  the  identity  operator^^^^^ 

In  the  first  part  of  this  dissertation,  we  develop  the  theo- 
retical foundations  of  our  approach  through  the  partitioning  technique 
extended  to  a multi-dimensional  reference  manifold. 

In  the  last  part,  we  apply  these  procedures  to  the  Stark 
effect  in  the  rigid  rotator  as  an  example  of  a one -dimensional  refer- 
ence manifold  and  to  the  two-electron  isoelectronic  series  as  an 
example  of  a multi-dimensional  reference  manifold. 


CHAPTER  II 


VARIATION  PRINCIPLE  AND  SEPARATION  THEOREM 


For  a Hamiltonian  operator  ~~3r\  which  is  bounded  from  below, 

a variational  calculation  gives  an  upper  bound  to  the  ground  state 

eigenvalue  of  'U.  The  variation  principle  states  cT = 0, 

A V 

where  ^ ^>Av  is  the  expectation  value  of  the  operator  : 


C-iJ?  > s y(dx) 

AV  “ jY*Y(dx) 


<yI  ^ Iy> 

< Y 1 Y> 


(2.1) 


Let  us  consider  that  the  trial  function  Y is  expanded  in  terms  of  a 
complete  set  of  functions  ) s (0^,02^ which  has  the  overlap 
matrix  with  the  elements 


A.  . 
ij 


= /0^*0j(dx) 


(2.2) 


The  Hamiltonian  is  then  represented  by  a matrix  H with  the  elements 


•Kij  = . (2.3) 


In  practical  calculations  the  complete  set  {0}  has  to  be  truncated  to 
a set  {0  = (0^,0^^  . . . ^0j^}  of  a finite  order,  say  order  n,  so  that 

n 


-I 

i=l 


C.0. 


1 1 


(2.1.) 


In  order  to  improve  our  expectation  value  using  this  finite  discrete 


3 


1+ 


set^  we  invoke  the  variation  theorem.  For  an  expectation  value  of  an 
operator  bounded  from  below^  this  variational  technique  gives  an  upper 
bound  to  the  eigenvalue  of  the  ground  state.  Upper  bounds  to  the 
excited  states  are  estimated  through  Rayleigh-Ritz  process  by  the 
expansion  of  the  trial  function  in  the  form  of  {2.k) . 

Denoting  the  column  vector  of  coefficient  C.  by  (E,  and  the 

t 

adjoint  of  C by  C , one  can  write  an  equivalent  equation  to  (2.l) 
in  terms  of  matrix  notation 


<^>AV 


= E. 


(2.5) 


Varying  the  coefficients  and  applying  the  variation  principle^  one 
obtains 


E -E^ij)  Cj  = 0^  (i=l,2,...,n)  . (2.6) 


The  roots  E^^  •••>  of  the  associated  secular  equation 


det  . - E A.  .)  = 0 

^ ij  ij 


(2.7) 


give  the  upper  bounds  to  the  n lowest  eigenvalues  of  in  order , 
If  we  introduce  a projection  operator  associated  with  the 


truncated  set  {0/}  , such  that 
n 


& <i)^  = 0^  , 


and  characterized  by  the  relations 


( ±~1,2, . . . ,n) > (2.8) 


0^1^  = 0^,  Cr^  = and  Tr(  &^)  = n, 


(2.9) 


y T 

■^The  symbol  f applied  as  a superscript  to  an  operator  or  to 
a matrix  denotes  its  Hermitian  adjoint. 
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then 


defined  by 


^n  ^ ^n^  ^n 


(2.10) 


gives  eigenvalues  •••}  ^ri<l  which  are  equal  to  the  roots  of 

the  secular  equation  (2.6).  Hence  (2.6)  can  be  replaced  by  an  eigen- 
value equation  of  the  form 


= E.  (5-V, 


111' 


(i+l,2,...,n)  , (2.11) 


with 


= U. 
1 1 


Since 


^ (^n^n  ^ 


(2.12) 


it  is  possible  to  diagonalize  the  JhH  and  ^ matrices  simultaneously 

by  a unitary  transformation  which  leaves  the  roots  of  (2.7)  unchanged. 

Therefore  any  theorem  or  property  related  with  the  roots  of  (2.j)  does 

not  depend  on  the  choice  of  the  basis  as  long  as  the  basis  spans  the 

same  space.  Let  us  further  introduce  an  additional  function  0 , ^ to 

n+1 

the  truncated  set  [0)^,  so  that  [0]  = [0^  ,0^, . . . ,0  ,0  , , } and  an 

n n+i  I d n n+i 

associated  Hermitian  projection  operator  Q-  characterized  by  the 
relations 

^n+l^n+1  ^ ^n+1^  Tr(  = n+1 , 

(2.13) 

^+l(^n  ^ ^n^n+1  ^ ^n  ^n+1 

n+1  n+1 

Then  defined  by 


-P?l 


n+1 


^n+1  ^ ^n+1 


(2.1k) 
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will  have  n+1  eigenvalues.  For  the  simplicity  of  notation,  we  will 
sometimes  use  m for  n+1  in  this  chapter. 

It  is  interesting  to  see  how  these  upper  bounds  given  by  the 
variation  principle  are  improved  by  enlarging  the  subspace  in  which 
the  trial  function  Y is  varied.  This  effect  is  qualitatively  explained 
by  the  well-known  "separation  theorem"^  . which  states  that  the  n 

eigenvalues  of  are  upper  bounds  to  the  n eigenvalues  of  in 
order  from  below,  and  the  eigenvalues  of  and  satisfy  the 

following  inequality,' 


E ^ h ^ E , T 
n n n+1 


(2.15) 


where  h.  denotes  the  n eigenvalues  of  and  E . denotes  the  n+1 

1 1 

eigenvalues  of  ;^n+l^ 

Let  us  now  study  the  separation  theorem  in  somewhat  more  detail 
in  connection  with  the  degeneracy  and  the  common  eigenvalues  of 
and  ^ since  they  are  closely  related  to  the  asymptotic  behavior 

and  the  lost  eigenvalues  in  the  partitioning  technique^^^ 

In  the  following,  for  simplicity,  we  will  work  with  the  ortho- 
normal set  of  {{5  } 


m 


m 


Let  us  consider  a matrix  |}-JJ  with  the  natrix  elements  of  the 


form 


<Afij  (i,j=l,2,...n,m)  ^ 

and  transform  the  matrix  by  means  of  the  unitary  transformation 


y 


(2.16) 
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where  (p  is  the  n-dimensional  unitary  matrix  consisting  of  eigenvectors 

n *2 

of  (HI  . Then  we  obtain  a secular  equation  of  the  form: 


where 


and 


<^mm 


V 


Im 


V, 


2m 


V 


nm 


V 


ml 


hl-E 


V 


m2 


hg-E 


i=l 


■mi  ik 


V = V* 
km  mk 


V 


mn 


= 0, 


h -E 
n 


\k  S 1^2,...,n), 


(2.17) 


Denoting  V = *V  , 

mk  mk  mk 

one  can  express  (2.I7)  by  the  form 


(K  -E) 


mm 


TT  (hi-E)  - f V^[  TT  (h.-E)] 
1=1  k=l  J 


= 0. 


(2.18) 


Defining  a continuous  function  Y(£.)  differentiable  to  any  order  by 

v(£).(«^-a  -|v2j|^{h.-e)]  , (2.19) 

(is'k) 

(2.18)  states  that  m eigenvalues  of  0^ are  the  m roots  of  Y(£:) : 


Y(E)  = (3^^-E).TT(\-E)-f;v^[T  (»i-E)]  ■ 0- 

1=1  k=l  1=1^ 

( i?^k) 


(2.20) 


*2 

Compare  W1 


th 


(8.20). 
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Let  us  assume  that  the  eigenvalues  of  are  arranged  in  non- 

decreasing order 


h,-  h^  ^ 


^ h . ^ h , 
n-1  n ’ 


then  the  behavior  of  the  function  Y(^)  is  conveniently  explained  by 
considering  the  various  possible  cases. 


Case  (l).  None  of  the  elements  V , are  vanishing  and  all  the  h.  are 
' — mk  ° 1 

different  from  one  another. 


In  this  case^  from  (2.20)  we  have 

■ - 4 ,4  C'l-h) 


(2.21) 


n-1 

Y(h  ) = - TT  (h.-h  ) 
' n'  mn  . , ^ 1 n 

1=1 


It  is  seen  from  (2.20)  that  Y(h^)  < 0,  furthermore  the  sign  of  Y(h^ 
is  different  from  that  of  Y(h^).  Hence,  there  exists  at  least  one  real 
root  between  Y(h^  and  Y(h^).  We  also  note  that 

Y(-“)  = lim  1^  - 1 ■fi’ 


(2.22) 

_)>  0,  if  m is  even 


^ 0,  if  m is  odd. 

Since  the  signs  of  Y(-“)  and  Y(h^)  are  different  from  each  other, 
there  exists  at  least  one  root  which  is  smaller  than  hj^  and  in  like 


manner  the  different  signs  of  Y(h^)  and  Y(-h»)  guarantee  that  there 
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exists  at  least  one  root  which  is  higher  than  h^.  Furthermore  (2.I9) 
cannot  have  more  than  m roots.  The  above  reasoning  gives  the  following 
separation  theorem. 

Theorem  2.1.  If  none  of  the  elements  of  V , are  zero  and  all  of  the  h. 
mk  1 

are  different  from  one  another,  (2.I9)  has  m distinct  roots  satisfying 
the  inequality, 

Et  < h,  < 4 h_  < < E < h ^ E . (2. 

1122  nnm  ^ 


Case  2.  All  of  are  nonvanishing  and  some  of  h^  are  degenerate. 

In  order  to  facilitate  our  discussion,  we  will  adopt  the  nota- 
Pi 

tion  h^  , where  the  superscript  p^  denotes  the  degree  of  degeneracy  in 
th 

the  i distinct  level  in  order  from  below,  so  that  the  following 
relations  are  satisfied: 


u ^ u 2 uPl  u ^2  . ^ u ^ 

h^  = h^  = . . . = h^  hg  ~ ^2  ^ “ “2 

= h^^  = ...  ^ (2.24) 

where 

h^  = h^,  h^  = hg,  . . . , h^  — hp^,  h^  - ^pj^+l  > (2.25) 

£ 

S Pi  = n . (2.26) 

. i=l 

Then  (2.I9)  will  have  the  form 


(2.28) 


10 


with 


j'  = 


j"  = 


k-1 

t 1 


i=l 

r p 


< ^ 

I Pi 

i=l 


and  (2.20)  takes  the  form 


j -E)J  (hJ-E)-  f [f  (hJ-E)J] 

11  1=1  k=l  /i=lx 

(i#k) 


for  k = 1, 
for  k > 2, 

for  k = 1^ 

for  k >,  2, 


(2.29) 


= 0.  (2.30) 


The  difference  between  (2.20)  and  (2.3O)  in  their  form  is  that  (2.3O) 
has  an  additional  factor  if  (h^-5)^Pi  and  hi  remains  as  a root  of 

i=l  1 1 

(2.27)  with  the  degeneracy  (p^-l).  We  can  use  the  same  reasoning 
applied  to  (2.I9)  to  Y(£) : 

Y(t) 


Y(£)  H 


(2.31) 


and  obtain  the  separation  theorem  of  the  form 


1 X , 1 - = p(Pl-l)  = uP  - , 1 _ „1 


^ h^  Eg  • • • 


hf  < E^  < . . . < h^  ... 

^2^  ^ ^2i+l 


(2.32) 


where 


^1  ^1^  ^2  ^2’  ^2 


Pl-1  1 

^ "Pl+1  = "3  ’ 


It  is  seen  from  (2.32)  that  the  change  from  ^ to  ^€. 
reduces  the  degeneracy  of  each  level  by  one  and 
erate  roots  between  the  degenerate  sets. 


m 


has  non-degen- 
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Case  (3).  Some  of  are  vanishing  and  none  of  the  h^j^  are  degenerate. 
If  an  element,  for  example,  is  zero,  we  can  write  (2.I9) 


and  (2.20)  as 


Y(£)  = (h.-z)[(-f{  -E)  T (h.-f)  - Z i Tf  (h..-£)3l 

- J " ^inin  i=l  k=l  nik  ( i=i  ^ x ^ 


(2.33) 


(i;^h  m) 


(2.3k) 


In  this  case  the  eigenvalue  h.  of  remains  as  an  eigenvalue  of  ^m^ 

n * 

This  is  due  to  the  fact  that  and  Z Uj^j  are  non-interacting  with 
respect  to  . Since 


Y(6)  ^ Y(£)/(h.-£:) 


(2.35) 


has  n distinct  roots  as  discussed  in  Case  (l),  Y( £ ) has  m distinct 
roots  if  the  j root  of  Y(£),  which  is  located  between  h^  ^.and 
is  different  from  h.,  and  the  following  inequalities  are  satisfied 
according  to  whether  the  j root  of  Y( 5 ) is  lower  or  higher  than  h^ : 


< h,  < ...  <E.  < h.  = E._^,  < h._^T  < E 

112  J J 


. , , ^ T X. . _ . . . h < E , (2.36) 

j+1  j+1  j+2  n m ' 


or 


E,  <:  h < E^  ...  <E.  = h.  4 E.,,  < h . _ <.  E . _ ...  <h  < E . (2.3?) 

112  J J J+1  J+1  J+2  n m v ^ 

t h 

If  the  j root  of  Y( £ ) happens  to  equal  to  h^,  then  one  degeneracy  of 
the  level  h^  is  added  and  gives  the  following  separation  theorem  : 

El  < hi  < E2  . . . Ej_l  C E^  = h^  = E^^^  < h^^^  < E^^^  . . . < h^  < E^. 

(2.38) 

The  above  discussion  can  be  extended  easily  to  the  case  where  more  than 


one  is  vanishing. 
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Corollary  2.1;  Degeneracy  may  be  induced  in  the  level  in  the  spec- 
trum of  only  when  the  corresponding  V . is  vanishing. 


Case  ( h) . Some  are  vanishing  and  the  corresponding  h^  is  degenerate. 

This  is  the  combined  case  of  Case  (2)  and  Case  (3)  and  (2.27) 

and  (2.33)  have  to  be  considered.  In  connection  with  (2.27),  we  can 

2 

consider  two  different  cases;  one,  is  vanishing^  and  the  other, 

2 2 

is  not  vanishing.  If  vanishes,  for  example^  one  can  write  (2.27) 
as 


Ji 


(i?^k,i?^j) 


(2.39) 


k^i 


and^  as  one  has  seen  in  the  discussion  of  (2.33),  degeneracy  for  the 
level  h^  will  remain  and  may  be  increased  by  one  through  the  change 
from  to  ■K";  and  if  is  not  vanishing^  then  one  can  go  back 

to  (2.27)  and  it  reduces  to  the  Case  (2). 

The  four  cases  discussed  above  give  a lucid  look  at  the  sepa- 
ration theorem  in  connection  with  the  degeneracy. 

Since  (2.39)  is  a special  form  of  (2.27)  (2.33).  in  ■ 

Chapter  III  we  will  not  consider  (2.39)  separately. 


CHAPTER  III 


BRACKETING  FUNCTION  AND  FIRST  ORDER 
ITERATION  OF  EIGENVALUES 


Using  some  of  the  equations  derived  in  Chapter  II,  we  will  now 
introduce  a function  f( £ ) called  bracketing  function.  Only  the  eigen- 
values of  different  from  those  of  are  bracketed  consecutively 

by  use  of  first  order  iteration. 

From  (2.20),  for  each  root  E of  Y(£'),  one  obtains 


(3.1a) 


for  the  case  that  all  roots  h.  are  distinct  and  all  V , are  non-van- 

1 mk 

ishing.  In  case  there  are  degeneracies  in  the  spectrum  of  ;^n 

one 

obtains  for  the  roots  of  Y( £ ) as  defined  in  (2.31): 


E = 


TT(h.-E) 
i=l  " 


(3.1b) 


and  for  the  case  V . vanishes,  one  has,  from  (2.35) 

m j ’ ’ \ i 


n 


E,  1 

k=l  ifk 

■ £ = (^^j) 

‘^mm  ^ n 

.E  (hi-E) 

(Si) 


(3.1c) 
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Ik 


We  note  here  that  (3* lb)  does  not  indicate  the  degree  of  degeneracy  of 
h^  and  (3-lc)  is  independent  of  h^(the  corresponding  is  vanishing). 

In  the  later  applications  of  (3-lb)  and  (3-lc)  to  the  iteration  pro- 
cess, where  0 is  the  reference  function  of  a one-dimensional  manifold 
one  has  to  note  that  there  is  no  information  about  the  degree  of  degen- 
eracy nor  about  the  eigenvalue  h^  corresponding  to  a vanishing  j • 

ro  1 

In  this  sense,  the  eigenvalue  of  h^  is  lost^  . With  these 
facts  in  mind,  we  will  primarily  work  with  (3-la).  One  can  write 
(3-la)  in  the  form 


E = f(  E,  ) 


(3-2) 


for  an  eigenvalue  E,  where'  the  function  f is  defined  by 

„ V2 


"a.- 


(3-3) 


for  the  continuous  variable  £.  Since  f'(  £.)  < 0 : 


S>f(S) 

c?£ 


V 


mk 


2l"  ('^k  ~ 


< 


(3.k) 


roi 

£.  and  f( £ ) will  bracket  at  least  one  eigenvalue  . In  this  sense, 

we  will  call  f( £ ) a "bracketing  function." 

It  is  interesting  to  derive  the  bracketing  function  by  use  of 
the  partitioning  process  and  we  will  introduce  L8wdin's^  ^ develop- 
ment of  the  partitioning  technique  in  terms  of  an  operator  formalism. 
Let  O'  he  a self-adjoint  projection  operator  which  defines  a certain 
subspace  S((?')  of  order  g in  the  total  Hilbert  space,  so  that 


15 


- e-\  Tr((9')=g*^  . (3.5) 

The  orthogonal  complement  to  the  subspace  S(  C?")  is  defined  by  the 
operator  P = 1 - <9^,  which  satisfies  the  relations 

p^  = p,  p'^  = p and  er?  ^ Yor  = 0.  '(3.6) 

Let  us  further  introduce  the  operator  T = T(£1)  by  the 
definition 

T(£)  = P[a.  (?'+P(£-  ^)P]'^P  ^ (3.7) 

where  o;  0 is  an  arbitrary  number.  In  the  following  we  will  often 
use  the  symbolic  notation 

^ ^ • ^3.8) 

It  is  seen  that  T fulfills  the  following  relations: 

= p . ‘ (3.9) 

e-^  = IS'  = 0 . (3.10) 


Let  us  now  consider  the  operator  -f\-  defined  by  the  relation 

_TL  ( £ ) = <^+ T(  d );C^  <9- . (3.11) 

Using  (3-9);  we  obtain  the  identity 

P(:^-<?)-A  = P(^ -£)(?"+  P(^-OT^e^=  «.P5£e^9'=0.  (3.12) 


*3 

In  practical  applications  of  the  partitioning  technique,  we 
often  utilize  the  bracketing  property  of  the  first  order  iteration  pro- 
cedure. In  these  problems  the  choice  of  an  arbitrary  subspace  S(<? ) 
has  some  importance  in  connection  with  the  convergence  property  of  the 
iteration  method.  In  this  sense  for  g = 1,  we  will  call  the  subspace 
S((9')  a one -dimensional  reference  manifold,  and  for  g^  2 , we  will 
call  it  a multi-dimensional  reference  manifold. 
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for  all  values  of  5 • Hence  we  have 

(7^-f)-n.=  (^+p)(^ -£)-a  = 0'{'^-B)-Cl. 


If_O.(E)0  is  an  eigenfunction  of with  an  eigenvalue  E,  0 has  to 
satisfy  the  relation 

-E)^(E)0  = 0 , (3.13) 


or 


E0^0  = . (3.U) 

This  is  the  necessary  and  sufficient  condition  for  E to  be  an  eigen- 
value of  . Let 


with 


g 


(T-  y I0.X0.I 

1^1 


Then  the  condition  (3.I4)  gives 

y {'ft . ,)  4 0 (0>  = 0 , (i=i,2,...,g) 

.4j  J 

J = 1 

where  "^^(E)  “ 5^  ^ 

^ij=<0j^l0.> 


(3-15) 


(3.16) 


(3-17) 
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From  (3.16)^  we  see  that  E has  to  satisfy  the  secular  equation 


^,1-E 

^12  • • ' 

■5^1 

. . . 

^2s 

^gl 

■^2  • • • 

A. 

(3.18) 


and  further  <5^0  cannot  be  arbitrary  if  Tr(<9")  = g > because  the  com- 
ponents < 0j|0>  must  be  eigenvectors  of  We  will  work  with 

(3.18)  in  more  detail  in  Chapter  V in  connection  with  the  multi-dimen- 
sional manifold. 

Let  us  now  define  an  operator  3k': 


■m:  - p^p  , 

and  its  eigenfunctions  so  that 

= h.^! 

^ X 11 

then 


:(E)  = E 


°°  I Y'  > < Y ' I 

i i 


(E  - h.) 
i=l  1 


( i 1,2, ... f) 


(3*19) 


,(3-20) 


(3-21) 


i=l  1 


(3-22) 


For  the  case  g = 1^  one  can  put  &=  with  1^  then 

(3.16)  reads  simply 


- E)<  f |0>  = 0. 


(3-23) 


18 


Assuming  < fU/i  > ^ 0,  one  obtains 


E = ^ ) cp> 


» <'?IKI?>  -h<<?  iKTfijFt  1?>  , 


The  introduction  of  (3 *21)  into  (3-24)  gives 

C e _ K.O 


E 


0<5» 

= <^|-B[q’>  -t-  Z 

Us.i 


(3.25) 


If  we  put 

^ mm  = <(  O’  I 1 ^ ^ 

v^.  . K |tl'>  ( 

then  it  is  seen  that  (3.24)  is  of  the  same  form  as  (3*3);  which  indi- 
cates that  the  function  f(£  ) defined  by 


f(£-)=  I -K  K>  4-  < 9 1 ^ Trs.;  ^ 19>  ^3.26) 


is  a bracketing  function  of  E. 

The  first  order  iteration^  using  this  bracketing  property^  will 
be  convergent  if  | ^ I ^ 1 ^ a-nd  will  be  divergent  if  j>  j 

In  this  connection,  the  choice  of  the  reference  function  is  important. 

In  the  treatment  of  a Hamiltonian  which  can  be  written  as 
the  sum  of  two  terms,  ^/  = + V,  it  is  convenient  to  introduce  a 

reduced  resolvent  T^  and  a reaction  operator  t through  the  relations: 


■ (£  - fifV 

t(£  ) - V + VTj^jV  . 


(3.28) 
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The  operators  and  t will  hence  depend  on  £ . Using  (3.26)  we  will 

try  to  find  a practical  method  for  obtaining  lower  bounds. 

By  using  the  relations  TV  = T t and  T = T + T tT  , LBwdin  has 

o o o o _ 

derived^"^^  the  following  transformation  for  the  operator 

(3.29) 


Substitution  of  this  expression  into  (3.26)  gives 

f(£)  = <<?  |«M^>  + ^ I ?> 

- I f<°l  <?>  + <9 1 1 1^> 


(3.30) 


where 


cf  = C I To  K ° ? 


In  connection  with  the  solvability  of  the  SchrBdinger  equation 
It  Is  convenient  to  take  as  a reference  function, 
since  for  this  case  (3-30)  simply  reads 

f(6)  = e:"  4.  <qp"’(  V/ -t  V \ (3-31) 

T>  . 

In  principle  we  can  obtain  the  upper  bound  £,^  to  any  desired 
accuracy  and  by  use  of  the  bracketing  property,  (3-30)  or  (3.31)^  we 
can  estimate  a lower  bound 


+ <9\ ±C£.> \<f>  , 

I <r">  ■*' <'7°l  T°>  . 


(3.32a) 


(3.32b) 
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However^  practically,  it  is  often  difficult  to  work  directly 
with  the  reaction  operator  t(£.),  and  we  are  obliged  to  develop  some 
method  which  gives  an  operator  t' {£.)  characterized  by  the  inequality 


. (3.33) 

Then  (3.32a)  and  (3.32b)  are  replaced  by 

■?“>  + <^i  {3-3i*) 

£[  E <r  I I , (3-35) 

and  we  have  the  following  inequalities: 

£",<d,<  E<  . (3-3b) 


A useful  way  to  achieve  this  end  is  to  introduce  the  concept  of  the 
"inner  projection" as  discussed  in  the  next  chapter. 

AS  B means  < Y)All^>  for  every  ^ in  the  common 

domain  of  A and  B. 


CHAPTER  IV 


THE  ONE -DIMENSIONAL  REFERENCE  MANIFOLD 
AND  LOWER  BOUNDS 

Ldwdin^^^^^  has  developed  a method  for  the  estimation  of  lower 
bounds  to  the  energy  eigenvalues  of  the  Hamiltonian  "^/^=  V.  In 

this  method,  in  order  to  avoid  the  ordering  theorem  of  the  "Intermediate 
Hamiltonian" ^ , the  bracketing  theorem  as  well  as  the  concept  of 
"inner  projection"  is  used. 

Let  A be  a positive  definite,  self-adjoint  operator.  Denoting 
the  eigenvalues  and  eigenfunctions  of  A by  Ct^  and  U.,  respectively, 
one  obtains  the  spectral  resolution 


A - 


(4.1) 


with  0.^  > 0. 


Let  us  further  define  an  operator  A^  by  the  relation 

A^  = E CL.2  I U.XU.l  ; 

i 

then  the  "inner  projection"  of  A,  A',  is  defined  by 


a'  = A2  q’^  A2  , 


(4.2) 


(4.3) 


where  is  a Hermitian  projection  operator  on  a linear  manifold  F^, 
which  is  spanned  by  a set  of  linearly  independent  functions, 

f H (f^,f2,...,f^).  (4.4) 
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Denoting  the  adjoint  of  £ by 

•£+  = 


r \ 


* 


n 


one  obtains  the  corresponding  metric  matrix 

= rt  = <f  If  > 

A typical  element  of  ZSi  is  defined  by 

Aw  = I f^>. 


(4.5) 


(4.6) 


(4.7) 


With  this  notation  one  can  define  a Hermitian  projection  operator  Q by 


-f  = l£> 


(4.8) 


with  the  property 


and  Tr(Q’")  = n . 


(4.9) 


Then  (4.3)  will  have  the  form 


A'  = A2  1 A2 


(4.10) 


If  we  incorporate  one  more  function,  say  into  the  set  and 

denote  the  corresponding  projection  operator  by  Q with  the  property 


+ 


Q»+l()"+l  . q“+1  . q"+1  and  Ir(Q“-"‘)  - n+1. 


(4.11) 
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[gl 

then  one  obtains  ^ 


Q^+1q^  = Q^Q^+1  = q’" 


(4.12) 


Equations  (4.9)^  (4.10)  and  (4.11)^  together  with  the  relations 


q“  = qV  = q"Vso  , 


(i-q“+44i-q“'"4  io  , 


(4.13) 


give  the  inequality 


0 &q"  i i X 


{4.14) 


or 


0 ^ A2q\2  £ A2q’^^^A2  ^ A 


(4.15) 


Thus^  by  enlarging  the  space  the  "inner  projection"  A'  approaches  A. 

In  the  following  it  is  convenient  to  introduce  three  more  mani- 
folds (g^.g^,  . . .,g^),  flq=  (h^,h2,...,h^)  and  j=  ( j jg, . . . , j^) 


defined  by  the  relations 

i:=A'^^,  f=A"^lK  and  f = A"'^(d-'R°)  J . (4.l6) 

Substitution  of  (4.l6)  into  (4.3)  gives 

. A'  = |A  (4.1?) 

with  1 A|  ; (4.18) 

A'  = (4.19) 

with  A = <hU~'^  lh>  ; (4.20) 


and 
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with 


(4.21) 

(4.22) 


The  "inner  projection"  in  the  form  of  (4.10)  is  called  by  LBwdin^^^  the 
"standard  inner  projection"  whereas  (4.17)  and  (4.19)  are  called  the 
"Aronszajn  projection"  and  the  "Bazley  projection,"  respectively. 

Since  the  "inner  projection"  of  the  form  (4.2l)  was  first  introduced 
by  Lhwdin^^^ , we  may  call  it  the  "Lhwdin  projection."  For  conven- 
ience we  will  also  call  the  manifolds  ^ , IK  a-^d  j the  "standard 
space,"  the  "Aronszajn  space,"  the  "Bazley  space"  and  the  "LBwdin  space," 
respectively  . 

According  to  (3.32a),  the  two  quantities  6^  and  £^,  which 
bracket  a true  eigenvalue  E are  connected  by  the  relation 


-^<^1  l^>, 

where  is  an  upper  bound  to  E and  is  a lower  bound  to  E.  Further- 
more if  one  can  evaluate  5 | which  is  smaller  than  £ 

^ ^ E < , (4.23) 

then  6 1 is  a lower  bound  to  E.  This  can  be  accomplished  by  introducing 
the  operator  t'  which  is  smaller  than  t.  If  t is  positive  definite  we 
can  always  obtain  t'  by  means  of  an  "inner  projection," 

t'  = t^Q\^  . (4.24) 


r?! 

Lowdin  has  shown  that,  even  if  t is  not  positive  definite,  we  can 
form  t'  such  that  t*  ^ t for  a positive  definite  perturbation  V. 
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Consider  an  energy  level  of  whose  upper  bound  satisfies 
the  inequality 


£ < E 


'p+1 


(4.25) 


where  is  the  (p+l)^^  eigenvalue  of  from  below.  Then  the 

choice  of  an  Aronszajn  space  for  the  inner  projection  of  V so  that 


is  sufficient  to  satisfy  the  inequality  t'(  £.^)  ^ t(  , 


where 


t'  = V’  + V'TV' 

= Iv-VT^  V)  V , 

if  ^ includes  all  the  lower  lying  unperturbed  eigenfunctions  Y?  with 


(4.26) 


respect  to  Hence  the  minimum  requirement  for  ^ is  given  by 


(4.27) 


This  choice  of  the  manifold  is  conveniently  applied  to  the  problems 
where  V/Jf  can' be  expressed-  in  terms  of  a known  finite  number  of  eigen- 
functions of>^^^^^. 

One  example  of  this  case  will  be  given  in  Chapter  VI. 

r 1 ry  1 Q1 

Bazley  and  Fox^  ^ have  discussed  the  similar  cases  in 
connection  with  their  special  choice  of  linear  manifold. 


CHAPTER  V 


THE  MULT I -DIMENS TONAL  REFERENCE  MANIFOLD 
AND  LOWER  BOUNDS 

The  main  difficulty  in  the  use  of  (3.34)  or  (3.35)  the 
calculation  of  lower  bounds  is  that,  for  an  upper  bound  greater 
than  the  first  excited  state  E°  of  the  generalized  reaction 

opertor  t(£^)  is  not  positive  definite.  Even  though  we  can  circum- 
vent this  difficulty  by  use  of  (4.26),  its  applicability  is  rather 
limited  due  to  the  appearance  of  the  reduced  resolvent  T^  in  t. 

These  two  difficulties  and  the  ordering  theorem  are  avoided  by  use  of 
the  multi-dimensional  reference  space  in  the  partitioning  technique. 

The  general  outline  of  this  partitioning  process  is  introduced  in 
Chapter  III. 

— 

Let  us  consider  the  operator  and  defined  by 

^(J)=  (5-1) 

and 

> <5-2) 

where  O'  is  a Hermitian  projection  operator  defined  in  (3.5)  with 
Tr  (gr)  = g.  Since  T(£)  is  a function  of  £,  so  is  and  -K,  and  g eigen- 
values  of'i^(£.)  will  change  continuously  with  £ except  at  some  singular 
points.  At  this  point,  it  is  convenient  to  consider  an  eigenvalue 
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problem  of  the  form 

/w' 

~ ■’  • • • )S)  t (5*3) 

It  is  seen  that  (3.114-)  is  a particular  case  of  (5 -3)  with  both  £ and 
6 £ equal  to  E which  is  an  eigenvalue  of  . 

Using  the  property  of  the  complementary  projection  operator  P; 

'P  = 1 - O' 


we  see  that 

= 0 > 

and  any  function  in  S(P)  *5  is  an  eigenvector  with  an  eigenvalue  zero 
for  all  values  of  6 ; this  case  is  of  little  interest  in  connection 
with  the  partitioning  technique^  and  we  will  confine  ourselves  to  S(&) 
so  that 


P0^  = 0. 


(5-^) 


Therefore  for  eigenfunctions  of  which  are  in  S(P),  this  partitioning 
process  fails  to  give  the  corresponding  eigenvalues. 


Theorem 


*1.1.  If  £ is  equal  to  one  of  the  eigenvalues  of 


i.e.,  6^(^)  then  is  an  eigenvalue  of  . 

Proof:  Starting  from  (5.3)^  one  obtains 


(5.5) 


^It  is  convenient  to  refer  to  the  subspace  associated  with  the 
P as  S(P). 
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o--H{^+  ^Ce,)-HO-)t  = ^^0-0 
^(K-  Tr^-)Kc^3  = o. 


Since 


Z(i^-^.)((^^  ="  by  (3.12), 

(6-  t)(H-^^)Co--\' 


or  (-K- C<5--b  ©-J)  - O j (5-7) 

this  indicates  that  ^ is  an  eigenfunction  of 

with  an  eigenvalue  6^.  Q.E.D. 

Furthermore  from  the  relations  (3*17)  ^o.d  (3*18)  there  follows: 
Corollary  5.1.  If  E is  an  eigenvalue  of  , but  not  an  eigenvalue 

r\^ 

of  P :f^p,  then  E is  also  an  eigenvalue  of 

In  Figure  I,  the  relation  between  the  eigenvalues  of  ( £) 
and  those  of  ^ are  schematically  illustrated  and  in  order  to  explain 
this  figure  we  will  discuss  the  asymptotic  behavior  of  the  eigenvalues 
of  £ ) as  £ varies  from  - «>  to  From  (3.22)  one  obtains  for 


(£)  = 

= 6t  ( R + ^ Tre;  ^ 3 


(5.8) 


where  ■u'=  P P 


and  = h^Y^ 
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^2^  ^3^  \ and  denote  the  eigenvalues  of  H; 

and  h denote  the  eigenvalues  of  Ml  • 
an  e^^  and  e^  denote  the  eigenvalues  of  rui  xhe 
author  is  indebted  to  Professor  P.  0.  Lbwdif for 
showing  this  diagram. 
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It  is  sometimes  more  convenient  to  consider  this  partitioning  process 
in  terms  o2  a secular  equation  of  a finite  dimension  and  later  extend 
this  to  the  infinite  Hilbert  space  to  match  with  the  operator  formalism 
shown  above. 

If  one  uses  an  orthonormal  basis  ,55^, . . . ) in 

g 2^  ■’  g 

expressions  (2.5)  (2.6) , one  obtains  the  matrix  equation 

IHl  C ^ EC  (5.9) 

with  , (5*10) 

and  Kl  . (5.11) 


Let  us  now  "partition"  the  basis  into  two  subsets  (a)  and  (b)  and 
denote  the  associated  subspaces  by  S(a)  and  S(b)^  respectively.  The 
matrix  ]H1  and  the  vector  (C  may  be  written  in  the  form 


1 iKlaa. 

IHOab  \ 

(C  = 

f u \ 

7 

H = 

and  instead  of  (5*9)^  two  equations  result: 

+IH2<x(j>Ci^=-  JrCa.  ^ 
C i,  -=■  H.  C b 


(5.12) 


(5.13) 

(5.14) 


If  the  eigenvalue  E of  IH]  is  different  from  the  eigenvalues  of 
(5.14)  gives 

= C ^ ' [HI  ba  C4  . (5.15) 
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Substitution  of  (5*15)  into  (5. 13)  leads  to  the  relation 


^ <13- a. 


(5.16) 


where 


■*"  IH)ab  ( ^ Xbb  IHJbb))  * IH^pci.  . 


(5.17) 


If  the  inverse  of  (E  J1  exists^  (5.I5)  gives  the  relation 

between  and  and  if  (5.I6)  can  be  solved,  can  be  obtained  by 
(5»15)»  However  for  an  eigenvalue  of  IH;  which  is  equal  to  one  of  the 
eigenvalues  of  THl^b^  this  partitioning  device  fails  to  give  the  eigen- 
functions as  well  as  the  eigenvalues,  and  in  this  sense  the  eigenvalue 
is  lost  (see  the  statement  after  (5*^))» 

Equation  has  exactly  the  same  form  as  (5*3)^  with  the 

understanding  that  S(b)  is  extended  to  the  Hilbert  space  of  infinite 
order  orthogonal  to  S(a).  Transforming  the  matrix  {J-Q  by  means  of 
the  unitary  transformation 


/ : 0 

\ 

\ o i y.. 

where  the  matrix  (ij^bb  composed  of  the  eigenvectors  of  IJ'^bb'’ 


(5.18) 


obtains  the  equivalent  matrix  IH)'  of  the  form  shown  in  (5*19) 


H'.  ! 

iHlb  \ 

djU  j 

J 

(5.19) 
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*6 


(5.20) 


where  are  eigenvalues  of  From  this  matrix  iO'D ' 

one  obtains  the  equivalent  form  of  (5.I7) 

pOi.x  = Maa  -t-  Mi,(  Elbb  - Mtb)  ‘ ^ (5.21) 

where  the  matrix  (E  the  diagonal  form  with  an 

element  7 — for  the  i^^  element.  Equation  (5.2l)  is  equivalent 
(E  - h^) 

to  (5.8). 

It  is  seen  from  (5.8)  and  (5.2l)  that  there  are  two  types  of 

.-v*.  — 

asymptotes  for  the  eigenvalues  of  -^/(£)(or  corresponding  to 

the  eigenvalues  of  <5-^^  ( or  J and  P;^P  (orlHQ^^),  the  former  is 
horizontal  and  the  latter  is  vertical.  The  number  of  asymptotes  of 
each  type  is  given  by  the  number  of  distinct  eigenvalues  of  il^JV’^and 
p -Uv,  respectively.  However  the  common  eigenvalues^  if  they  exists 
of  and  P 3^  P are  eliminated  here. 

In  order  to  examine  the  behavior  of  the  eigenvalues  of  ^ (£) 
(or  H (£))  qualitatively,  an  example  is  considered  where  the  dimen- 
sion  of  S(a)  is  two  and  the  dimension  of  S(b)  is  three.  The  operator 
formalism  and  the  matrix  formalism  will  be  used  interchangeably 
according  to  convenience.  Consider  the  non-degenerate  case  of 

^We  assume  here  that  the  unitary  transformation(5 . 18)  is  so 
arranged  that  h^  ^ ^2  ~ ^3  ’ * * * 
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and  denote  the  eigenvalues  of  by  e^  and  e (with  e^  •<.  e ) and 

diSi  L ^ 1 2^ 

those  of  by  h^^  h^  and  h^^  respectively. 

Assuming  S(a)  is  spanned  by  0^  and  0^  S(b)  is  spanned  by 
0^^  0|^  and  0^,  (5.21)  has  the  form 

Paa(S^  = Maa  + ^ 4b  ‘ Mb)'"  ML  (5-22) 


with 


H 


aa 


and 


-^12  \ 

V ^21  ^2/ 

/ K.l 

' 

11 

t - 

.t-  = 3 

\ S-  Hh 

e- 

Note  that  h^  in  Fig.  1 is  related  to  by 

1 JJ 

" ”^33’  ^2  " '^W  ^3  " "^55  • 


(5.23) 


(5.24) 


Defining 


Hla(£)  " Kb(tAb  - Kb>'  > 

(5.22)  reads 


(5.25) 


Maa(^)  = H.  ^ «:,(£)  . 


(5.26) 


In  what  follows^  the  horizontal  asymptotic  behavior  of  the 
eigenvalues  of  is  discussed  in  Case  (l)  and  Case  (2)^  whereas 

SL3l 

the  vertical  behavior  is  studied  in  Case  (3). 


34 


For  the  horizontal  behavior  of  the  eigenvalues  of  ( £)  one  has 

3i3. 

C ± ^ +- ce;  = (5.27) 

where  Is  a matrix  each  element  of  which  goes  to  zero  as  ^ . 

It  is  well  known^^^-’^^^^^  that  if  operators  A and  B bounded 
from  below  have  a common  domain  and  satisfy  the  inequality 

A 'C  B , (5*28) 


then  the  eigenvalues  of  A are  smaller  than  the  eigenvalues  of  B in 
order  from  below.  Furthermore,  we  can  denote  the  inequality  of  the 


matrices  A and  IB: 

A < IB 

if  every  diagonal  element  A^^  of  A 

diagonal  element  of  B^^  of  tB  , i.e., 

A.  . < B.  . 
11  11 


(5.29) 

is  less  than  the  corresponding 

(5.30) 


for  every  arbitrary  choice  of  basis  in  that  domain. 

The  character  of  I}"13  (£•)  as  a function  of  £ is  determined  by 

aa(£)  • Since  the  numerator  of  every  diagonal  element  is  greater 
than  or  equal  to  zero,  for  an  arbitrary  basis  chosen  in  S(a),  the  sign 
of  each  diagonal  element  is  determined  by  the  sign  of  corresponding 
denominator. 


Case  (1).  - <C£<h^ 


For  this  case  every  diagonal  element  of  (£)  is  negative 
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for  an  arbitrary  choice  of  basis  in  S(a)^  and  the  inequality 


M (i)  < (5.31) 

is  satisfied.  Furthermore  every  diagonal  element  of  for  an 

arbitrary  choice  of  basis  in  S(a)  becomes  smaller  as  ^ increases  its 
value  in  this  range.  Consequently^  the  eigenvalues  of  IK!(£.)  are 
decreasing  monotonically  from  their  corresponding  horizontal  asymptotes 
e^  and  e^  as  one  increases  the  value  of  £ . In  this  sense  the  hori- 
zontal asymptotes  are  upper  bounds  to  the  eigenvalues  of  (£) 

3.3. 

for  £ £ h^ . 


Case  (2)  . h^  ■^  £ < + 00 
5 

For  this  range  of  <£  , every  diagonal  element  of  (£)  for  an 

3.3. 

arbitrary  choice  of  basis  in  S(a)  is  positive^  giving  rise  to. the 
inequality 


> THl 


aa 


(5.32) 


and  the  eigenvalues  of  IHl  (£)  are  decreasingly  monotonically  and 

3.3. 

approach  their  corresponding  horizontal  asymptotes  e^  and  e^,  respec- 
tively. Hence  the  horizontal  asymptotes  are  lower  bounds  to  the  eigen- 
values of  ]!HJ  (£)  for  £ > h . 

3.3.  5 

Case  (3)’  £ approaches  h^(£-S'h^). 

In  order  to  facilitate  consideration  of  this  limiting  case,  the 
operator  formalism  expressed  in  (5*8)  will  be  used: 


(5.33) 
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where 


= f 

i=l 


•(£)  ^ (£-hJ 


For  generality^  the  dimension  of  S(b)  is  not  restricted  but  in  the 
present  example^  Tr(P)  = 3.  Rearranging  (5*8)  one  obtains 


eiiit)  »-‘d- 


'<<p.  o 

(i^k) 


(£- V 


Defining  the  operator 

(£-\) 

and  denoting  its  eigenfunction  and  eigenvalue  by  ^ and 

plkS-  5k  2 

one  obtains 


(T. 

(5.3^) 

(5-35) 

(5*36) 


2 

k* 


= 


A-lk 


Combination  or  (5 *37)  ^o.d  (5-38)  gives 


5k  [5k 


(£  - \) 


0, 


I 


or 


k = 0 , 

<<1^ 


^k  = 


(£  - hj 


(5.37) 


(5-38) 


(5.39) 
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Therefore  S,  has  two  distinct  eigenvalues,  zero  and  K I ^ ^ 1 

C £ - 

/ _ 

k “ ■ 

If  ^ k ~ using  the  relations  C9y  = 0^  and 

PY'  = Y'  one  obtains  : 
k k-’ 

= [ Kte]  Etfe 

- C«- 

Hence  is  a simultaneous  eigenfunction  to  and  P^P  with  the 

eigenvalue  h^.  Since  the  inverse  of  P(<£  - :;V)P  does  not  exist  for 

hj^^  the  partitioning  technique  fails.  Furthermore  h^  is  no  longer 
a singular  point  of  (5 -34)  and  (5 '35);  accordingly  there  is  no  asymp- 
totic behavior  at  £•  = h^^  and  the  eigenvalue  is  lost  in  this  sense. 
Therefore  it  will  be  assumed  that  f 0.  Under  this  assumption 

there  exists  at  least  one  non-zero  eigenvalue  The  eigenvalue 

whether  it  is  zero  or  non-zero^  may  be  degenerate^  however^  one  can 
always  diagonalize  in  S(^). 


unless  O 


In  order  to  facilitate  the  following  consideration^  let  the 

subspace  S(  (^)  be  divided^^^^  into  orthogonal  subspaces  S((3'^)  and 

Q. 

Ic  \ I -fc 

S((^,J  according  to  the  two  distinct  eigenvalues  of  such  that 


C£  “ 


(5.40) 
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and 


H le  " H k 


= Ot  P = 0 , 


(5.41) 


where 


0^+  0-b 


0-^0-  = <?■  = 


ol  ot  = o'^  o::  = o 

C?'V>  0~ 


ot  . (5-'*2) 


Introducing  the  relations  (5*40,  (5-41)  and  (5-42)  into  (5*34)^ 
one  obtains 


= 0-i:  i^r  S iJ  o-^  ^ 0-^<^ 


(5-43) 


(l 


where 


€fo  H fe  + 5,  c £-1,.) 


(5-44) 


It  is  convenient  to  introduce  the  operators  jsf(^)  and  i^(<£),  ^ 

3.  D 


by  the  definitions^  as  in  (3*7)  ' 
k r-  ^ ^ ^ k 


H . (5.45) 


JCOi,  3 (5. 

where  /3  is  an  arbitrary  non-zero  number. 

In  the  following  we  will  simply  use  the  symbolic  notation 


46) 


Jcc)^  - 


C (5 


(5.4?) 
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C s - ^c<L>  1 


(5.48) 


In  order  to  observe  the  asymptotic  behavior  of  & as  £.  approaches 
hj^^  one  can  contract  with  respect  to  the  subspace  S(^^)  and 

S(<9'^)^  respectively,  using  (5>43)  and  (5.44),  one  obtains  (cf.(3.l7)) 


d-t^,  ol  = 

(5.49) 

* 

(5.50) 

where 

= 6L  K O'  + 

(5.51) 

(5.52) 

As  £ approaches  h^,  ^ 
operator  (£.)  go  to 

3 


becomes  enormously  large,  which  makes  the 
k 

zero.  Hence  one  has  the  results 


(5-53) 


(5-54) 


£-5’  4^, 


These  give  the  theorem  that,  if  £.  approaches  one  of  the  eigenvalues  of 

say  h^,  where  h^  is  not  an  eigenvalue  of  the  n^  eigenvalues, 

where  n^  = Tr(  of.  Q ^ approach  the  eigenvalues  of  G(<r)^^ 


and  the  n eigenvalues,  where  n = Tr((?^^),  of  <?<^^become  infinite 
^ 3.  3. 

(T  £- 


with  the  order  of  magnitude 


4o 


In  concluding  the  discussion  in  connection  with  Figure  1,  the 
way  of  estimating  lower  bounds  to  the  ground  state  will  be  studied  in 
connection  with  (3- 18). 


Theorem  5.2.  If  is  the  ground  state  energy  of  and  the  eigen- 
function of  the  ground  state  is  not  orthogonal  to  S(  0~) , then  the 


lowest  eigenvalue  of 


is 


Proof:  By  corollary  5*1^  has  at  least  one  eigenvalue 

which  is  equal  to  Ej^.  If  one  denotes  the  eigenvalues  of  ^(E^)  by 
Ef,  E^,  in  order  and  assumes  that 


E 


1 


k > 1, 


(5-55) 


then  by  the  variation  principle  one  has 

^1  ^ ^ 

where  h^^  is  the  ground  state  eigenvalue  of  E'ff'2.  Since  the  cases 
where  the  eigenvalues  of  are  equal  to  the  eigenvalues  of  P^P, 
are  excluded  from  consideration,  one  may  assume  that  E^^  ^ hj^,  then 

by  (5.55) 


E,  h. 


(5.56) 


Now  consider 


C E«  - (5.5-j.) 

The  term  2j  — i T (5*57)  is  negative  definite  by  the 
i=l  i 

inequality  in  (5.56),  therefore  if  one  considers  the  range  Ej^  1^^^^^ 


4l 


then 

= -^{\)  < ^(£)  , (5.58) 

and  it  follows ^ that 

'e . ^ <f^  (i=l,2,..^  Tv{(^))  , (5.59) 

th 

where  ^ is  the  i eigenvalues  of  ^^(£,).  While  decreasing  £.  con- 

-«W 

tinually  from  Ej^  to  each  will  increase  in  magnitude  continuously, 

and  it  will  finally  reach  the  point  where  which  is  greater  than 

is  equal  to  6.  Then  by  theorem  5.I,  the  £,  at  this  point  ( £ = 

is  an  eigenvalue  of  , however  E^^  = E^^  and  this  contradicts  the 
assumption  that  Ej^  = E^^  if  k > 1.  Q.E.D. 

One  may  assume  that  the  Hamiltonian  ^ can  be  written  in  the 

form 

■A’-  #f“  + v 

where  ^o 

has  known  eigenvalues  and  eigenfunctions  and  V is  positive 
definite.  It  is  further  assumed  that  ^o 

has  discrete  energy  levels 


E°  ^ E°  £ ...  , (5.60) 

below  its  continuous  spectrum.  The  corresponding  orthonormalized 
eigenfunctions  are  denoted  by  T?,  so  that 


' ' 1 1 i 


< T?(  T?  > 
■ J 


ij 


(5.61) 


If  the  projection  operator  in  (3*5)  is  associated  with  the  functions 

...o  , 

by 


g 

= S fy?  X Y°  I 
i=l  ^ i 


(5.62) 
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the  introduction  of  (5.6I)  into  (3.l4)  leads  to  the  form 


^ c|>  *=  C5-  C -h  V + VTce^;  V ■)  O'  £^. 


(5.63) 


It  is  convenient  to  introduce  the  generalized  reduced  resolvent  and 
the  reaction  operator  t for  the  multi-dimensional  reference  manifold 
by  the  relations 


(5.64) 


t(e.)  s V + VT(£)V 


(5.65) 


Then  using  the  operator  identity 


(a-b)"^  = A"^  + A'^B(A-B)"^  , 


it  can  be  shown^^^  that 


t-1  = V-'  - T 


(5.66) 


(5.67) 


provided  V exists. 

The  condition  which  makes  T^(£)  negative  definite  will  be  the 
sufficient  condition  for  t ^(6.)  and  t(S)  to  be  positive  definite.  If 
the  variable  S is  below  the  continuous  spectrum  of  with  the 

relation 


(5.68) 


then  T^(£)  is  negative  definite  if  one  uses  the  projection  operator 
defined  by 


i=l 


(5.69) 
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where  ^ i k.  This  gives 

t > 0 , (5.70) 

and  by  means  of  the  inner  projection  defined  in  (4.3)  one  obtains 

t'  = t2  q"  , (5-71) 

so  that 

t'(£)  i t(£)  . (5.72) 

In  order  to  consider  the  bracketing  property  for  an  eigenvalue  E of  the 
Hamiltonian  in  a multi -dimensional  reference  manifold,  consider 
(3.14)  and  define  a function  f(£)  by 


then 


^ (5*73) 


<?')«>> 

d€ 


(5.74) 


It  is  seen  from  (5.74)  that  for  a properly  chosen  reference 

r i2i 

function  0,  which  satifies  (3.I4),  the  bracketing  theorem  holds^ 

hence  there  will  be  at  least  one  eigenvalue  E of^  between  € and  f(E)j 

therefore  if  is  an  upper  bound  to  the  eigenvalue  E and  if  both  E 

*7 

and  c- ^ are  between  two  consecutive  distinct  eigenvalues  of  P , 
then  f(£  ) is  a lower  bound  to  E: 

^1=  • (5-75) 

The  first  order  iteration  procedure  using  this  bracketing  property  will 
be  convergent  if  and  divergent  if  > !• 

**^Note  that  the  common  eigenvalues  of and  P'^  P,  if  they 
exist,  are  eliminated. 


Assume  that  denotes  the  ground  state  eigenfunction  of 
with  the  eigenvalue  and  further  assume  that  an  upper  bound  to 
is  lower  than  the  lowest  eigenvalue  of  P which  is  not  a simul- 
taneous eigenvalue  of  •ff’.  If  satisfies  the  relation 

o 

p o p+1  ' 

then  the  introduction  of  a projection  operator  of  the  form  (cf.  (5»69)) 


^ - Z |y°XY°I  , (where  g p) 

i=l  ^ ^ 


(5.76) 


makes  the  reaction  operator  t(6^)  positive  definite.  Hence  all  eigen- 
functions of  associated  with  the  eigenvalues  less  than  the  upper 

bound  are  included  in  the  set  { Let  us 
o 1 2 p g' 

define  by  the  relation 

> (5-77) 

where  (S'  is  given  by  (5.76).  Using  the  definition  (5.77)  and  relations 
(5-63).  (5.65)  and  (5.73)  f one  can  write 

£(4„)  - . (5.78) 

For  this  case  the  relation  in  (5*75)  ™ay  be  replaced  by 

^l<  E^<  (5.79) 

Defining  by 

= 4 0^1  ^°  + t'(£.^)l0j^>  /<  0^1  85^7,  (5.80) 


and  noting  the  relations  (5*72)  and  (5«79)  one  obtains  the  following 
inequality : 


(5.81) 
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Since  is  not  known,  one  minimizes  the  operator  , 

-U"  - -ff°+  f(.tj  , (5.82) 

in  S((5^  to  obtain  some  value  less  than  <S|  . 

Denoting  the  lowest  eigenvalue  of  in  S{&-)  by  one 


obtains 


Then  it  is  seen  that  is  given  by  the  lowest  eigenvalue  of  the 
Hermitian  matrix 

\ 


(5.83) 


hi 

‘la 

...  tj 

Ig 

f 

^0  f 

t 

21 

^2  *^22 

• • * *^2g 

• • • • 

1 

f 

^0 

gi 

*'g2 

...  E + 1 

g 

(5.84) 


/ 


Combining  (5*79)  (5*83)  one  obtains 

£U  /-  ' 

1 -^1  i '=l<  • 


(5.85) 


This  method  will  be  applied  in  Chapter  VII  in  the  estimation  of  the 
lower  bound  to  the  ground  state  energy  of  two-electron  isoelectronic 


series . 


CHAPTER  VI 


THE  STARK  EFFECT  IN  THE  RIGID  ROTATOR 

In  this  chapter,  lower  bounds  to  the  energy  eigenvalues  for 
the  rigid  rotator  in  a uniform  electric  field  (the  Stark  effect)  are 
calculated  using  the  method  described  in  Chapter  IV.  For  the  inner 
projection  of  the  positive  perturbation  V,  the  Aronszajn  space  ^ is 
used,  where  ^consists  of  the  eigenfunctions  of  and  has  to 

satisfy  the  condition  given  by  (k.2j) . The  notation  is  consistent 
with  that  of  Chapter  IV.  The  equation  used  in  this  chapter  for  the 
calculation  of  lower  bounds  is  (3.3^+)  with  t'  expressed  by  (k.26): 

?>  + <^lvl^x^lv  - VT^Vl^>"^<^lVlf >.(6.1) 

In  a problem  where  V is  a positive  perturbation  with  a con- 
stant coefficient  K as  a strength  parameter,  one  can  write 

V = KU  . (6.2) 

For  this  case  one  obtains 

e;  - . (6.3) 

where 

B = <^|U/K  - UT^U|^>  , (6.4) 

<f\v]  9 > 
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(6.5) 


If  is  an  eigenfunction  of  , then 

¥ = ? 


(6.6) 


The  wave  equation  for  the  rigid  rotator  consisting  of  two  mass 
points  and  which  are  separated  by  fixed  distance  R,  is 


2 

s in^,  e 5 

where  W is  the  rotational  energy  and  M is  the  reduced  mass  of  mj^  and 


tiH  h ® h’>  * TT-  + w°(9,0)  = 0, 


r ^ 

2MR^ 

[sin  0 

dd 

1 . e • 


M = (m^m2)/(mj^+m2) 


(6.8) 


Let 


-(u(0,0)  = 


iir?  h k'>  * 


sin  9 


(6.9) 


and 


(2MR^(jo)/fi^  = E° 


(6,10) 


where  o)  is  the  Hamiltonian  in  dimensionless  units  and  E*^  is  the  corre- 
sponding energy.  If  denotes  the  normalized  wave  function  given  by 


[20] 


o ^ r2j+l  (^-  Iml  .)!„  , q) 

U+lml):  ■'  ^ 


im0 


(6.11) 


where  P^(cos  ©)  is  the  associated  Legendre  polynomial,  then  (6.7)  reads 


simply 


k8 


or 


(6.12) 


The  rigid  rotator  of  dipole  moment  ji  in  a uniform  electric 
field  F is  characterized  by  the  wave  equation  of  the  form 

= , (6.13) 

with  "pf  = CJD  + K.cos  6 . 


and  K = (2uFMR^)/ft2  . (6.1U) 

In  order  that  the  perturbation  V in  (6.2)  shall  be  positive  definite, 
we  shall  write 


(JO  - K 


(6.15) 


V = K(1  + cos  0)  , (6.16) 

and 

U = 1 + cos  0 . (6.17) 

Then  this  division  of  ^ into  an  unperturbed  and  a perturbed  part 
gives 

+ KU  (6.18) 

with 


U >0  . 


(6.19) 
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For  this  case,  the  SchrcJdinger  equation  reads 


[cu(0,0)  + K cos  0]i'(e,0)  = . (6.20) 

This  equation  is  easily  seen  to  be  separable  with  respect  to  the  vari- 
ables 0 and  0,  so  we  can  write  the  eigenfunction  in  the  product  form 


T(e,0)  = z(e).  §(0) 


(6.21) 


with 


•>6 


i(«)  - ( 4;^) 


where  m = 0,  1,  2,  . . . ; 


and  Z(e)  = H C (cos  e)  , (6.22) 

i=0 

where  is  an  associated  Legendre  polynomial. 

Since  the  space  with  which  we  are  concerned  can  be  subdivided 
according  to  the  value  of  m,  we  can  treat  the  problem  separately  for 
each  subspace.  The  normalized  reference  function  0 in  (6.1)  is  arbi- 
trary; however,  in  order  to  make  the  first  order  iteration  of  (6.1) 
convergent,  the  condition^^^ 


<T#e9)TRf>  < <q>|<p>  (6.23) 

must  be  satisfied.  The  left  side  of  this  inequality  approaches  zero 
as  K approaches  zero  for  a fixed  reference  function  only  if  0 is  an 
eigenfunction  of  pc  ; therefore  it  seems  reasonable  to  make  this 

*8 

Except  for  the  case  m - 0,  we  have  doubly  degenerate  eigen- 
states, so  we  will  consider  only  the  absolute  value  of  m. 
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choice  for  0 provided  K is  small.  Thus  in  the  subspace  m,  we  have  for 
,th 


the  £ excited  state, 

. N>”(cos  0) 


(6.2h) 


where  Y^(e,0)  is  the  normalized  "spherical  harmonics"  and  where  is 
the  normalization  factor  for  the  spherical  harmonics  Y^(9,0)  given  by 


m r (^-m):(2^+l) 


[ 


(6.25) 


(^+m)I  liir  J 

Lower  bounds  for  the  energy  levels  can  be  calculated  by  use 

of  (6.3)  through  (6.5)  with  £,  in  T obtained  from  a variation  calcu- 

o o 

lation.  The  Aronszajn  space  which  satisfies  the  condition  of  {h.2j) 
is  given  by  , 


with 


, = y"* 

’i  nrf-i 


(6.26) 


Since  U = 1 + cos  0, 


^ = <gilu(Y|J>  = <y;J^._^((1+cos  e)fY|J>.  (6.27) 


Using  the  recurrence  relation 


(cos  0)P^  = (2m)‘^[(m-rn)p“^^  + (>g-Hn)p™_^]  , (6.28) 


one  obtains 


P.  = 

1 


(^+l-m)(^+l+m) 


(2m)  (2^+3) 


'^nrfi-l,m  ^ ®nrfi-l,i 


(^+m)(^-m) ■ 

(2m)(2<g-l) 


I 

®m+i-l,i-l  ; 


therefore 
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^ r (^+1-1 

L (2^+1 


sKj+l+ra) 

i)(2i+3) 


Pi  = 1 


P.  = 
1 


(^+m)(^-m) 


(2i+r)(2l-iy 


P.  = 0 

1 


for  i = 


for  i = i+l-m, 


for  i = H-m, 


otherwise, 


(6.29) 


and  = < Y^J(co-K)Iy^>  = ^(m)  - K . 


(6.30) 


we 


In  order  to  have  at  least  one  non-zero  element  for  P.,  as 
can  see  from  (6.29),  n must  be  at  least  as  large  as  (i-m);  otherwise  . 

£ [ = Z(m)  - K , (6.31) 

which  is  trivial. 

For  the  calculation  of  matrix  elements  of  (p,  we  consider 

■ <8il(U/K)  - UT^Ulgj>  . (6.32) 

h ■ (6.33) 


(1  + cos  S)C3.1  , 


where 


P = 1 - |<f  ><9/=  1 - Iy“><y^1 


Using  (6.28),  we  obtain 


i+1 


^ = D y"* 

1 , K nrfK-1 


with 


(6.34) 


(6.35) 


D, 


. , = 1 r (i-l)(i-2m+l)  12 

1-1  [£q-  (i+m-2)(i+m-l)+K]  L (2i-^2m-l!K2i+2m-3)J 

= 0 (if  .e  = i+m-2) 


(6.36) 
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D.  = — 

^ I ^o~  (i+m-l)(  i+m)+K] 


- 0,  (if  £ = i-hn-l) 


i+1 


[ ^o“(  ^"^)(  i+nrt-l)  ] 

0, 


if i+2m) 


(2i+2m-l)(2i+2nrfl) 

( if  i = i+m) , 


(6.38) 


Note  that  the  reason  vanishes  for  certain  values  of  i is  that  is 
orthogonal  to  Y^. 

Introducing  (6.33)  into  (6.32), 

Bki  = <gkl(U/K))g^>  -<gj^lu)F^>  . (6.39) 

Then,  the  non-vanishing  matrix  elements  of  IB  are 

B = f ~1  1 r i(  i+1 ) ( i+2m)  ( i+2nrH ) 1 2^ 

i+2,i  LT£  - (i+m)(i+m+l)-HCi(i+2nrfl)J  L (2i+2m-l)(2i+2nrf3)  J ’ 


i+l,i 


l-K  ISn"  ( i+m)(  i+m+l)+K^  fcn- 


■Q-  (i+m)(i+m+l)+K^  (Cg-  ( i+m-l)  ( i+m)+K 


t-('2i 


i(  i+2m)' 


(2i+2m+l)(2i+2m-l) 


]] 


B = i _r  if  i+^)  J r 1 ] 

K L (i-hn)(i+m+l)+Ki(2i+2m-l)(2i+2m+l)J  ■L£,^-(i-hn-l(i+m)-H<J 

- (izl)(i-l+2m) 1 

LICq-  (i+m-2)(i+m-l)-tK}(2i+2m-l)(2i+2m-3)  j 

1-1,1  [k  ((f^.  ( i+m-l )(i+m)+K)  (i+m-2)(  i+m-l )+k}] 

r f i-l)fi+2m-l)  I 2 ^ 

L(2i+2m-l)(2i+2m-3) J 

B = r ll- 1 rf  i-l)f  i-2)f  i+2m-l)f  i+2m-2)l2 

i-2,i  L {fg-  (i+m-2)(  i+m-l )+K}(2i+2m-3)JL  (2i+2m-l)(2i+2m-5)  1 * 

(6.1tO) 
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Equations  (6.29),  (6.3O)  and  (6.U0)  give  us  the  necessary 
matrix  elements  for  calculation  of  lower  bounds  expressed  by  (6#3). 
It  is  easily  seen  that 


®i+2,i  ®i,i+2 


®i+l,i  “ ®i,i+l 


If  we  change  the  Hamiltonian  in  (6.I3)  to 
^ = o)  “ K cos  0 


(6.41) 


we  would  expect  the  same  result  since  we  only  reversed  the  uniform 
electric  field.  The  algebraic  identity  of  these  two  cases  has  been 
demonstrated  as  a partial  check  on  the  algebra. 


Upper  and  Lower  Bounds  for  the  Rigid -Rotator--Numerical  Results 

The  entire  calculation  was  done  on  the  IBM  709  at  the  Univer- 
sity of  Florida  Computing  Center.  Upper  Bound  energies  were  first 

obtained  by  the  Rayleigh-Ritz  variational  method  using  subroutine 
■Xn 

Givens  (single  precision).  In  order  to  clarify  the  cases  where  the 

upper  bound  and  lower  bound  were  so  close  that  it  was  difficult  to  say 

which  was  lower,  the  upper  bounds  were  refined  using  a double-precision 
• • ^10 

Iteration  method  and  the  remainder  of  the  calculation  was  also 
carried  out  in  double  precision. 

*9 

"Eigenvalues  and  Eigenvectors  by  the  Givens  Method,"  Quantum 
Chemistry  Program  Exchange,  QCPE  12C  (I963). 

■^0 

C.  E.  Reid,  "Eigenvalues  by  Lbwdin's  Partitioning  Method," 
QCPE  l4C  (1963)^  modified  for  double  precision. 
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Except  for  the  level  m = 0,  all  the  other  levels  are  degenerate. 

However,  the  Hamiltonian  commutes  with  L = t therefore  one 

z 1 o (O’ 

can  separately  solve  the  energy  eigenvalue  problem  for  each  subspace 
for  a given  value  m of  L^.  Calculations  for  the  first  six  subspaces 
(S^  to  S^)  have  been  carried  out  for  the  twenty  lowest  eigenvalues  in 
each  subspace  with  the  perturbation  coefficient  K ranging  from  0.1  to 

1.0. 


In  Table  1,  upper  bounds  E = £ and  lower  bounds  E = £•'  of 

U o Jj  JL 

energy  eigenvalues  are  given  to  sixteen  significant  digits  for  typical 
choices  of  K and  m.  The  eigenvalues  are  labeled  by  the  Z of  the  un- 
perturbed state. 


Ny  is  the  number  of  basis  functions  used  for  finding  upper 


bounds  (these  functions  are  the  normalized  spherical  harmonics 

m’ 

^^1^1’  -1^’  number  of  basis  functions  ^ , 


^m+N  1^  Aronszajn  space, 

L 


This  method  of  evaluating  lower  bounds,  using  the  bracketing 
property  of  (6.1),  requires  good  upper  bounds,  and  the  upper  bounds  are 
improved  by  increasing  the  number  of  basis  functions  Ny  used  in  the 
variational  calculation.  The  results  in  Table  1 indicate  that  fairly 
good  upper  bounds  are  generally  obtained  by  taking  Ny  = + 5.  In 

order  to  obtain  good  lower  bounds,  we  have  to  enlarge  the  Aronszajn 
space,  and  we  can  see  from  the  results  indicated  in  Table  1 that  in 
the  present  case  the  dimension  of  the  Aronszajn  space  N_  is  sufficient 
to  give  good  lower  bounds  provided  N = i + 5. 

Li 


a 

0 

1 

2 

3 

h 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 
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TABLE  1 

UPPER  AND  LOWER  BOUNDS  FOR  EIGENVALUES  FOR 
THE  STARK  EFFECT  IN  A RIGID  ROTATOR 


m = 0, K = 0.^ 

\ ^U  \ " 


5 

5 

-0.410507934822997 

-0.410507934823054 

6 

6 

2.024369028254599 

2.024369028254599 

7 

7 

6.005966607238925 

6.005966607238923 

8 

8 

12.00277848867649 

12.00277848867648 

9 

9 

20.00016234906589 

20.00016234906587 

10 

10 

30.00106840547832 

30.00106840547831 

11 

11 

42.00075758568145 

42.00075758568144 

12 

12 

56.000565611185280 

56.00056561485278 

13 

13 

72.00043859831220 

72.00043859831218 

14 

14 

90.000035011^09683 

90.000035011(09682 

15 

15 

110.0002860416819 

110.0002860416819 

16 

16 

132.0002380955195 

132.0002380955193 

17 

17 

156.0002012884138 

156.0002012884136 

18 

18 

182.0001724138988 

182.0001724138987 

19 

19 

210.0001493429597 

210.0001493429596 

20 

20 

240.0001306165556 

240.0001306165555 

20 

20 

272.0001152074046 

272.0001152074045 

TABLE  l--Continued 


m = 0, 

K = 0.5 

n 

E and  E ® 

U Li 

17 

20 

20 

306.0001023751246 

306.0001023751242 

18 

20 

20 

23 

342.0000916098325  ^ 

342.0000915749897 

342.0000915751046 

19 

20 

20 

24 

380.0016458519702 

380.0000722250626 

380.0000822706509 

m = 1, 

K = 1.0 

1 

5 

5 

1-950333905229181 

1-950333905229141 

2 

6 

6 

6.011577927242672 

6.011577927242670 

3 

7 

7 

12.00832614435840 

12.00832614435839 

8 

8 

20.00551933162358 

20.00551933162356 

5 

9 

9 

30.00384626739139 

30.00384626739137 

6 

10 

10 

42.00281392300602 

42.00281392300600 

7 

11 

11 

56.00214127777589 

56.00214127777586 

8 

12 

12 

72.00168130591093 

72.00168130591091 

9 

13 

13 

90.00135388552912 

90.00135388552910 

10 

14 

14 

110.0011129663826 

110.0011129663825 

11 

15 

15 

132.0009307395922 

132.0009307395921 

12 

16 

16 

156.0007896715423 

156.0007896715422 

13 

17 

17 

182.0006782886897 

182.0006782886897 

TABLE  1- -Continued 


m = 1, 

K = 1.0 

E,,  and  E^  ^ 
U ij 

14 

18 

18 

210.0005888386502 

210.0005888386500 

15 

19 

19 

240.0005159358681 

240.0005159358679 

16 

20 

20 

272.0004557472937 

272.0004557472681 

17 

20 

20 

306.0004054860207 

306.0004054860206 

18 

20 

20 

23 

342.0003630874514 

342.0003630874365 

342.0003630874366 

19 

20 

20 

23 

380.0003274705419 

380.0003269927098 

380.0003269958156 

20 

20 

20 

24 

420.0062377999262 

420.0002179805039 

420.0002942543069 

m = 3 

II 

0 

-4i 

3 

5 

5 

11.99319648972468 

11.99319648972467 

4 

6 

6 

19.99888486422010 

19.99888456422010 

5 

7 

7 

30.00020901290087 

30.00020901290086 

6 

8 

8 

42.00053019645272 

42.00053019645271 

7 

9 

9 

56.00057406214888 

56.00057436214887 

8 

10 

10 

72.00053726950879 

72.00053726950875 

9 

11 

11 

90.000118038806258 

90.00048803806256 

10 

12 

12 

110.0004230273525 

110.0004230273523 

TABLE  l--Continued 


m - 3, 

K « 0.7 

£ 

\ 

E and  E ® 

U J-i 

11 

13 

13 

132.0003712115186 

132.0003712115184 

12 

14 

14 

156.0003262415718 

156.0003262415716 

13 

15 

15 

182.0002877983403 

182.0002877983402 

14 

16 

16 

210.0002550776513 

210.0002550776511 

15 

17 

17 

240.0002272074346 

240.0002272074345 

16 

18 

18 

272.0002033918627 

272.0002033918626 

17 

19 

19 

306.0001829503520 

306.0001829503519 

18 

20 

20 

342.0001653171585 

342.0001653172584 

19 

20 

20 

380.0001500277721 

380.0001500277720 

20 

20 

20 

420.0001367024596 

420.0001367024584 

21 

20 

20 

462.0001251239667 

462.0001250304234 

22 

20 

20 

506.0027336495416 

506.0000914202348 

m = 5, 

K = 1.0 

5 5 

6 6 


7 7 

8 8 

9 9 


29-99359063883720 

5 29.99359063883718 

14-1.99761861361275 

6 I41. 99761861361275 

55.99923213730367 

7 55.99923213730366 

71.99992678740871 

8 71.999926787^869 

90.00023337531917 
90.00023337531912 
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TABLE  1- -Continued 


m = 5, 

K = 1.0 

E and  E ^ 

U J_i 

10 

10 

10 

110.0003640283063 

110.0003640283062 

11 

11 

11 

132.0004112436340 

132.0004112436339 

12 

12 

12 

156.0004180542290 

156.0004180542288 

13 

13 

13 

182.0004054534825 

182.0004054534824 

14 

14 

14 

210.0003840228910 

210.00038140228908 

15 

15 

15 

240.0003591944730 

240.0003591944728 

16 

16 

16 

272.0003337620168 

272.0003337620166 

17 

17 

17 

306.0003091323675 

306.0003091323674 

18 

18 

18 

342.0002859711714 

342.0002859711713 

19 

19 

19 

380.0002645455813 

380.0002645455812 

20 

20 

20 

420.0002449100770 

420.0002449100769 

21 

20 

20 

462.0002270087367 

462.0002270087364 

22 

20 

20 

506.0002107319706 

506.0002107319655 

23 

20 

20 

552.0001962023146 

552.00019594^3010 

24 

20 

20 

600.0049841862676 

600.0001226802245 

lower 

^he  upper 
bound . 

bound  is 

given  first,  followed  by  the 

In  those  cases  where  more  than  one  N was  used, 
a lower  bound  is  indicated  for  each  N . 

Li 
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Upper  and  lower  bounds  agree  to  fourteen  significant  digits*^ ^ in 
those  cases  where  = ^ + 5.  For  the  highest  energy  levels 

listed,  for  example,  K = 1.0,  m = 1,  i = 20,  N = 20,  N = 20,  the 

U Li 

agreement  is  poorer,  but  is  improved  somewhat  by  increasing  N . The 

j_i 

limitation  here  seems  to  be  due  to  the  poorness  of  the  upper  bound 
rather  than  the  dimension  of  the  Aronszajn  space. 

To  see  in  more  detail  the  effect  of  N for  a fixed  value  of  N , 
several  examples  are  given  in  Table  2.  It  is  seen  from  (6.29)  that, 
when  the  Aronszajn  space  ^ is  chosen  as  indicated  in  (6.26),  for 
less  than  (<^-m),  the  lower  bound  Ey  is  a persistent  lower  bound, 
given  by 

Ey  = - K . 

In  order  to  obtain  better  lower  bounds,  therefore,  the  Aronszajn 
space  ^ has  to  be  chosen  in  such  a way  that  there  are  non-vanishing 
elements  of  ; namely  ^^-m+2  case.  Increasing 

Ny  it'  Table  2 beyond  the  values  indicated  did  not  improve  the  lower 
bound. 

*11 

Except  for  one  case,  m = 0,  ^ = 0,  K = 0.5,  where  it  agrees 
only  for  I3  significant  figures. 
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TABLE  2 


EFFECT  OF  ON  IN  THE  CONVERGENCE 

TO  THE  EIGENVALUE 


m 

K H 

Nl 

E and  E 

U Li 

0 

0.1  10 

15 

^ a 

110.0000114416483  . 

1-9 

109 . 9000000000000X-' 

10 

109.9477443614383 

13 

110.0000114416482 

0 

0.5  10 

15 

1-9 

110.0002860416819 
109 . 5000000000000^ 

10 

109.7412282838218 

13 

110.0002860416817 

15 

110.0002860416818 

1 

0.7  16 

20 

1-14 

272.0002233160454 
271 . 3000000000000* 

15 

271.6219989276613 

18 

272.0002233160452 

2 

1.0  12 

15 

1-9 

156.0007432191651 
155 ■ 0000000000000* 

10 

155.4101468398408 

13 

156.0007432191644 

15 

156.0007432191651 

3 

0.7  10 

12 

1-6 

110.0004230273525 
109 . 3000000000000* 

8 

110.0003141668453 

10 

110.0004230273523 

5 

1.0  20 

20 

1-14 

420 . 0002449 100770 
419 . 0000000000000* 

15 

419.3717170996393 

18 

1120.0002449100768 

20 

420.0002449100769 

from 

a_. 

The  notation  I-9 
one  through  nine. 

indicates 

that  the 

value  of  Nj^  ranges 

b. 

Asterisks  indicate  persistent  lower 

bound . 

CHAPTER  VII 


LOWER  BOUNDS  TO  THE  GROUND  STATE  ENERGY 
LEVELS  OF  TWO-ELECTRON  IONS 


The  non-relativistic  Hamiltonian  of  this  system  is  divided 
into  two  parts: 


+ 


V 


(7.1) 


with 


and 


T,  1 . . . 

V = — - in  atomic  units, 

h2 


(7-2) 


where  r^^  and  r^  denote  the  radial  positions  of  electrons  1 and  2;  r^^ 

denotes  the  distance  between  electron  1 and  2;  and  is  the  Laplacian 
operator  in  the  coordinates  ^^^(1  = 1^2).  Since  we  estimate  only  the 
ground  state  which  is  singlet  S,  we  restrict  ourselves  to  the  subspace 
in  which  the  spacial  coordinates  r^  and  are  symmetric  with  respect 

F21  1 

to  their  interchange,  and  the  angular  momentum  is  zero^  . 

In  the  application  of  the  method  developed  in  Chapter  V,  we 
have  freedom  in  choosing  a subspace  for  the  inner  projection  of  the 
reaction  operator  t(£).  However,  we  have  to  take  into  account  two  fac- 
tors in  our  consideration:  one  is  the  convenience  of  our  treatment  and 

the  other  is  the  convergence  property  of  our  lower  bounds.  Here  we 

introduce  two  types  of  spaces  for  our  inner  projection:  one  is  the 
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practically  convenient  Bazley  space  the  other  is  the  Lbwdin  space 
Since  the  choice  of  the  Bazley  space  has  some  interesting  features  in 
connection  with  the  simplicity  of  our  calculations  and  the  independence 
of  our  lower  bound  to  the  upper  bound  we  will  briefly  describe  this 
choice  in  Appendix  II,  while  in  this  chapter  we  will  mainly  work  with 
the  LBwdin  space. 

The  LBwdin  projection  of  t gives 

f il(vJ<°)  (7.3) 

where,  using  (4.22)  and  (5*67)^  A has  the  form 

A = < J I ( ^°)  - (V 

(7-4) 

The  unperturbed  Hamiltonian  fV  has  known  eigenvalues  and  eigenfunc- 
tions which  are  products  of  hydrogen-like  wave  functions.  The  discrete 
eigenvalues  in  atomic  units  are  given  by 


~ 2"  = 1,2,...) 


(7.5) 


Arranging  these  in  a non-decreasing  sequence 


. . . , 


(7.6) 


one  obtains  from  (7*5) 


E°  = 
^1 


E°  = - 0.625Z^,  E°  = - ^2, 


E°  = - O.53I25Z  , E°  = - 0.52Z^ 


(7.7) 


Hence,  if  an  upper  bound  to  the  ground  state  energy  lies  below 
2 

-O.625Z  , then  the  one -dimensional  reference  manifold  can  be  used. 

In  Table  3 it  is  seen  that,  except  for  the  case  Z = 1,  we  can 
use  (3*35)  the  lower  bound  calculations  of  two-electron  ions.  The 
following  six  eigenfunctions  of  are  used  for  the  multi-dimen- 

sional reference  manifold: 

'^'6  = (§t^32(^l)^2l(^2^^2(^°®  ®12^^  ^ 


where  are  hydrogenic  normalized  radial  wave  functions. 

[2^1 

Six  terms  of  the  Hylleraas  series  are  used  for  the  Lbwdin 
space.  The  normalized  forms  of  these  six  terms  are: 


4 = (2r),^,  -ns 


8n 


4 = (2^),^. 

'2  l6n^/6 


j ^2 


9611 


■T]S 


i =-^ 


s e 


-T]S 


(7.9) 


4 = j2nP  2 -ns 
h 96n  " ® 


(2n  5 2 -PS 

2f-  = ^ e ' 

° l60it  v3 


where 


s = r^  + r^,  t = r^-r^  and  u = r^^  ’ 


(7.10) 


One  can  enlarge  the  Lbwdin  space  in  order  to  improve  the  lower  bounds. 


TABLE 
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The  arbitrary  scaling  parameter  t]  is  varied  to  maximize  the  lower 
bounds  obtainable. 

It  is  seen  from  (7.7)  and  Table  3 that  for  T ( £ ) to  be  nega- 

0 0 

tive  definite  S(  has  to  contain  at  least  four  lower  lying  eigen- 
functions of  , namely  and  for  Z = 1,  whereas  only 

one  eigenfunction  is  sufficient  for  Z 2.  2.  If  one  uses,  as  a typical 
example,  all  six  functions  in  (7.8)  for  S(O'),  then  from  (7.3)  and 
(7.i+)  one  obtains 


i=l 

(7.11) 


and 


where 


^ ^ ' (7*12) 


Substituting  (7. 12)  into  (5*84),  then  the  lowest  eigenvalue  of  (5.84) 
is  a lower  bound  to  the  ground  state  of^.  Once  the  A.  matrix  is 
obtained  the  rest  of  the  calculation  is  elementary.  Most  of  the  overlap 
integrals  of  the  type  complicated  in  their  algebraic  form. 

Some  typical  examples  of  the  terms  which  comprise  each  element  of  the  ^ 
matrix  are 


= 3;28125(£q+ti^)^  . ^ 

+ 8.75(£o-Hi^)(Z-ti)  + (4  ^n  2 - 4.5)t1^(Z-ti) 

+ 12^^  2 - 29)t)3  + (6  ^n  2 + h.^)r\{Z-r])^  , 

<32  J2  > " ^ i*5ii(z-n)  , 


67 


Numerical  Results  and  Discussion 


The  numerical  calculations  was  carried  out  on  the  IBM  709 

electronic  computer  at  the  University  of  Florida  Computing  Center. 

In  the  estimation  of  the  lower  bound  to  the  ground  state  energy,  by 

■X-Q 

use  of  (5*83)^  Givens  subroutine  was  used. 


The  upper  bounds,  denoted  by  are  taken  from  Table  3 and 
the  lower  bounds,  denoted  by  are  optimized  with  respect  to  the 
scaling  parameter  tj.  It  can  be  seen  that  the  scaling  of  the  exponential 
coefficient  for  the  reference  function  does  not  introduce  any  new 
feature  in  our  problem.  In  the  following  tables,  N denotes  the  dimen- 
sion  of  the  reference  manifold  associated  with  the  eigenfunctions 

’ whereas  Np  denotes  the  dimension  of  the  Lbwdin  space, 

R 

^2^  * ’ ■ ^ ^Np  ’ energy  is  expressed  in  atomic  units. 

In  Table  the  lower  bounds  to  the  ground  states  of  two  elec- 


tron ions  from  Z — 1 to  Z - 10  are  listed.  The  upper  and  lower  bounds 
agree  up  to  two  significant  digits  except  for  Z = 1. 


TABLE  4 

UPPER  AND  LCfWER  BOUNDS  TO  THE  GROUND  STATE  OF  TWO- ELECTRON  IONS 
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In  Table  5^  the  effects  of  and  on  the  lower  bounds  for 
each  Z are  demonstrated  with  optimized  t]  at  = 6 and  = 6.  It  is 

seen  form  Table  5 that  the  lower  bounds  do  not  improve^  as  one  might 
hope  from  the  discussion  of  Chapter  V,  when  one  increases  the  value  of 
This  might  indicate  that  even  if  the  lower  bound  in  (5.78) 
should  be  improved  by  increasing  N , the  corresponding  improvement  in 
S'^  may  be  offset  to  some  extent  by  the  tendency  for  S.”  to  be  depressed 
as  the  size  of  its  secular  equation  increases  (according  to  the  varia- 
tion principle).  On  the  other  hand  the  increased  value  of  Np  gives  rise 
to  an  improvement  of  the  lower  bounds.  For  Z = 1^  the  most  outstanding 
change  comes  as  we  increase  Np  from  3 to  which  corresponds  to  adding 
to  the  LBwdin  space  spanned  by  and  whereas  for  Z 2 more 

than  50  percent  of  total  improvement,  from  Np=  1 to  Np=  ^ takes  place  at 
the  change  of  Np  from  2 to  3*  Hence  the  choice  of  the  proper  function, 
rather  than  increasing  the  number  of  functions,  for  the  LBwdin  space 
seems  to  be  more  significant  in  our  problem.  The  effect  of  the  scaling 
parameter  q on  the  lower  bounds  for  Z = 1 (with  Np  = 4,  5,  and  6)  and 
Z = 2 (with  Np  = 1,2, ...6)  is  illustrated  in  Figs.  2 and  3.  In  both 
figures,  the  increase  of  Np  up  to  6 does  not  improve  lower  bounds  appre- 
ciably. For  some  part  of  the  curves  the  improvement  of  the  lower  bounds 
by  the  increase  of  Np  is  too  small  to  be  seen  in  the  present  figures. 

As  a numerical  demonstration,  for  Z = 1,  q = 0.7  one  has 


4 

6 

4 

-0.5561906 

-0.5561905 

5 

-0.5560375 

-0 .5560366 

6 

-0.5560233 

-0.5560211 
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TABLE  5 (Continued) 
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Fig.  2.  Lower  bounds  (Ej^)  to  the  ground  state  of  h"  as  a function  of  . 

scaling  parameter  (t])  with  Nj^=4.  The  number  attached  to  each 
curve  denotes  the  value  of  N . 
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Fig.  3.  Lower  bound  to  the  ground  state  of  helium  as  a function 
of  scaling  parameter  with  N =1.  The  number  associated 
to  each  curve  denotes  the  value  of  N„. 


Ih 


For  H , the  lower  bounds  assume  a flat  top  type  curve.  The  flat  top 
behavior  of  the  curves  for  Z = 1 may  explain  the  large  gap  between 
the  upper  and  the  lower  bounds.  This  behavior  may  indicate  that  the 
functions  we  have  chosen  for  the  LBwdin  space^  although  suitable  for 
helium,  are  probably  not  entirely  suitable  for  h"  with  respect  to  the 
variation  of  t]. 

A good  choice  of  the  LBwdin  space  is  essential  in  order  to 
improve  our  lower  bounds.  In  contrast  with  the  behavior  for  Z = 1, 

^^8*  3)  Z = 2,  shows  a fairly  good  convergence  of  the  bound  as  we 
increase  Npj  furthermore,  the  curves  are  of  an  approximately  parabolic 
appearance  indicating  that  the  choice  of  the  LBwdin  space  is  satisfac- 
tory with  respect  to  the  variation  of  the  scaling  parameter  t]  . Interr 
estingly  enough.  Fig.  3 is  very  similar  to  that  obtained  by  J.  G.  Gay^^^] 
He  obtained  the  lower  bound  to  the  ground  state  of  helium  by  solving  the 
determinant ; 

det  -0.  (7.13) 

The  comparison  is  made  in  Table  6 between  the  lower  bounds  obtained  from 
r fc  1 

(T-Ib)  Gay  and  those  obtained  from  (5.84).  It  is  seen  from 
Table  6 that  the  lower  bounds  obtained  as  eigenvalues  of  the  intermediate 
Hamiltonian  are  slightly  better  than  those  obtained  by  the  present 
method,  but  the  differences  in  the  lower  bounds  of  those  two  methods 
decrease  as  one  increases  the  value  of  until  at  Np  = 6 there  is 
essentially  no  difference.  The  dependence  of  the  lower  bounds  on  the 
scaling  parameter  t]  for  Z = 4,  6,  8 and  10  is  shown  in  Fig.  4,  where  N 

R 

and  Np  are  fixed  at  the  value  of  6.  As  in  the  case  of  Z = 2,  all  curves 
are  roughly  parabolic. 


COMPARISON  OF  THE  LOWER  BOUND  ENERGIES  OF  PRESENT  METHOD  FOR  HELIUM  WITH  THOSE 
OBTAINED  FROM  EIGENVALUE  PROBLEM  OF  THE  INTERMEDIATE  HAMILTONIAN 
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Fig.  4.  The  lower  bound  to  the  ground  state  of  two  electron  ions 
for  Z = h,  6,  8,  and  10  as  a function  of  t]  at  N =6  and  N. 
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The  upper  bound  of  the  H is  improved^  by  O.OOII5  a.u., 
from  -0. 52531  to  -0.52646^  by  extending  the  Hylleraas-type  wave  func- 
tion from  three  terms  to  six  terms,  whereas  is  seen  from  Table  5 the 
corresponding  improvement  for  the  lower  bound  is  about  O.II3.  This 

shows  that  the  lower  bounds  improve  more  rapidly  than  the  upper  bounds 

*12 

do  by  introducing  more  terms  . Hence  one  can  hope  to  improve  the 
lower  bounds  of  the  ground  state  of  H by  introducing  carefully  chosen 
functions  in  the  Lbwdin  space. 

In  concluding  the  discussion  it  can  be  noted  that,  as  one  can 
see  from  Fig.  2,  the  most  outstanding  improvement  of  the  lower  bound  to 
the  ground  state  of  H takes  place  when  one  incorporates  the  function 
^1).  1^1°  l4e  Lbwdin  space.  The  introduction  of  the  function  also 
shows  remarkable  improvement  on  the  lower  bound  even  if  it  does  not 
render  significant  contribution  between  the  range  t]  = 0.68  and  T)  = 0.®4. 
In  terms  of  hydrogenic  radial  functions,  and  assume  the  form 

^4  ~ ^1^^21^ ^2^  ^ ^ ^ (7*l4) 

h ~ ^5 *-^10^ ’^l ^^32^ ^2^  ^ 

(7.15) 

where  and  are  normalization  constants,  showing  that  and  are 
independent  of  62,2^  angle  between  electron  1 and  2.  Since  R^^(r), 


om  one 


^21^^^^  ^ad  ^1^2^^^  have  different  special  distributions  fr 
another,  it  seems  that  the  introduction  of  some  functions,  which  can 
give  a more  flexible  radial  distribution  with  respect  to  the  variation 


*12 


See  Table  1 of  [I5] . 
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of  over  the  wide  range  of  space  ^ can  improve  the  lower  bound  satis- 
factorily. In  this  connection,  it  might  be  wothwhile  to  incorporate 
into  the  LBwdin  space  some  carefully  chosen  radial  functions  which  have 
two  different  orbital  exponents,  for  example,  a function  of  the  form 
investigated  by  Chandrasekar^^"^^ , 


- ( , -bri-aro  v/  , , 

-(e  ^ ^ + e ^)(i  + c*rj^2 


+ 


(7.16) 


It  may  be  also  worthwhile  trying  to  use  Hylleraas  series  given  in  (7.9) 
by  giving  a different  scaling  parameter  for  each  function. 

For  practical  applications,  a proper  choice  of  the  functions 
with  suitable  scaling  parameters  for  the  inner  projection  is  to  be 
preferred  over  expanding  the  number  of  functions. 


APPENDIX  I 


Evaluation  of  the  Matrix  Elements  used  in  Chapter  VII 


For  S-state  configurations,  the  atomic  wave  function  can  be 
described  in  terms  of  the  three  coordinates^^^^ . 


s = ^1  + r^  , 


- ri  - , 


(Al.l) 


u = r 


12 


The  volume  element  in  terms  of  s,  t,  u coordinates  (Hylleraas  coordi- 


nates) 


IS 


[28, p. 1737] 


dVfdv2  = fr^u(s^-t^)  du  ds  dt 
with  the  range 


(A1.2) 


U^t£.U,  Oiu:£s 


(A1.3) 


In  terms  of  the  Hylleraas  coordinates,  given  by  (7-2)  assumes 

the  form^^^^ 

= - 


2 >2  ,2 


■i;  + ^ ^ ^ gs(u"-A  , 2t(s^-u^) 

3s  8t2  <5u  u(s^-t^)  <3sdu  u(s^-t^) 


2 2.  ^2 


9tdu 


+ 


Us  d 


_Ut d_  + ^ A.  + UZs 


and 


(s^-t2)  3s  (s2-t2)  dt  u du  (s2_t2)j  ^ (A1.4) 

(AI.5) 


V = i =i 

ri2  u 


The  multi-dimensional  reference  functions  in  (7.8  ) are  linear  combi- 
nations of  the  form 


- _ -as  -bt  k ^ m 
f - e e s t u , 


(A1.6) 
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where  a>lbl>0  and  I,  is  even;  whereas  the  function  for  the 


inner 


projection  have  the  form 


J - e St  u . 


(A1.7) 


Matrix  elements  involving  functions  of  type  f and  j were  all  evaluated 
in  a closed  form.  In  the  evaluation  of  A matrix  using  Hylleraas 
coordinates  s,  t and  u,  two  types  of  integrals  are  utilized.  One  is 


(o:/n,k,^,m)  = /“  e‘<^s*^ds  /®u®du  /^(s^-t^)"^  t^-^dt,  (A1.8) 


where  n -1,0,1,  k ^0,  m ^-1  and  q;>0.  The  algebraic  expressions 

for  (A1.8)  are  given  in^^*^^.  The  other  is 


00 

(a,b/k,i,ra)  - e"^®s’^ds  /Qu"*du  £^(s^-t^)  t'^e”^'"  dt 


/l-I  y (q+k);  q 

.^+m+2^^p!  q- 

b p=0  q=0  ^ 


ifq+k+2)(q+k+l) 


(a+b)^'*’^'''^  ( (a+b)^ 


_ (it2)(i+l)1  (-1)"^'*  f(q+k-^g)(q+k■n)  Xii£iLe±Ull 

b"  J (a-b)2  ■ Jj 

£+2  / , ^,p+m 


■ . f-hn+b  -^1  2 

b p-^+1  ^ q=0  ql  l(a+b) 

• • {'  * '-"1 


^^+m+2  k+1 
b a 


U+2)iJl+l) 

b^ 


(k+2)(k+l)l  f (p-hn)I  , (iWl):(2^-hn+4) 

a^  -'Ip^O  P'  J £[ 


(A1.9) 


The  integral  (AI.9)  was  used  as  a means  of  checking  some  algebraic  expres- 
sions for  the  matrix  which  were  obtained  by  use  of  the  polar  coordi- 


nates 


[31] 
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For  the  convenience  of  Fortran  programming  this  form  is  transformed 
into 


p=0  qt'o 


) ( q+k+g ) ( q+k+1 ) 


(a+b)‘''"''+l  I (a+b)2 


Lmlimll  + (-1)°*^  /fq+k-t^Hq+k+l)  U.+2)n+l)'< 

I (a-b)'‘+>'+l  I (a-b)2  J 


r. 


, ^+m+2  _ 

b q=nrf2 


y (q+k)!b^  A (p+m)I 

.iUo  q'  pi 


p=q-m 


,(a+b) 


q+k+1 


f(q+k+2)( q+k+1)  _ (i-l^)(^+l)\  , C-l)"^q  /(q+i,+g)(q+t,+i) 

[ (a+b)^  b^  J 


(a-b)^+^+^  1.  (a-b)2 


U-^g)(^+l)^ 


-\  C^+m+l)  (’2^+m+4)  (q+k)!b‘^r  1 

7J’  b^^^  qib  ‘I-  l<a+b)'l^^‘^^ 


(a-b)'l'^*^+\, 

+ (l  + (-1)'"^ 


(^+k+n+2)! 
, 2 


i-iy 


_(a+b)^'^^'^3  (a-b)^’^^+3 


kl 

, ^+m+2  k+1 
b a 


^ ^ (p+m)  i f(.Ji+2)U+l) 


p=0 


I 


_ U-hn+l)ir2i;-hiri-U) 


(Al.lO) 


In  order  to  avoid  the  effect  of  round-off  error^  (AI.9)  is  programmed 
in  double  precision.  This  Fortran  Program  is  listed  on  the  following 


pages . 
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FUNCnUN  DHH{  A,B,K  ,L,M) 

C FUNCTEON  DCr(M)  AND  DCTD(M,N)  NEEDED 

C L,  M,  AND  N MUST  BE  ZERO  OR  POSITIVE  INTEGER. 

c B cannot  be  zero,  a is  positive,  the  absolute  value  of 

B IS  LESS  THAN  A. 

K1  = K+  I 
K2  = K f 2 
LI  = L + 1 
L2  = L + 2 
Ml  = M 4-  1 
M2  = M + 2 
LK  = L M 

KM  = K ♦ M 

LMl  ^ L + M 1 
LM2  = L + M + 2 
LM^  = L + M + 4 
KMl  = K + M + 1 

KM2  = K + M *■  2 

KLM2  =K+L+M+2 
KLM3  = K*L  + M*-3 
L2M4  = 2*L  + M + 4 
0 FK  = K 

D FL  = L 

D FM  = M 

0 FKl  = K1 

D FK2  = K2 

D FLl  ^ LI 

D FL2  = L2 

0 FKMl  = KMl 

D FKM2  = KM2 

0 FL2M4  = L2M4 

D ASQ  = A*A 

D BSQ  = B*8 

D AI  = l./A 

D El  = l./B 

D Also  = l./ASQ 

D BISQ  = l./BSO 

D APB  = A 4-  8 

D ANB  = A - B 

D BPA  = l./APB 

0 BNA  = l./ANB 

D APBSO  = APB*APB 

D ANDSQ  = ANB*ANB 

D BPASO  = BPA*BPA 

0 BNASQ  = BNA*BNA 

D FLLB  = ( (FL+2. )*(FL+1. ) )/  BSO 

D FKKA  = ( (FK+2. )*(FK+l. ) )/ASQ 

5 IF(L-{L/2)*2)  6,  6,  8 
D 6 SIGL  = 1. 

7 GO  TO  9 
D 8 SIGL  = -1. 

9 IF(  M-IM/2)*2  ) 10,  10,  12 
0 10  SIGM  = 1. 

11  GO  TO  13 
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D 12  SIGM  = -1. 

0 13  RCl  = DCr(L)/(B**LM2) 

0 14  Q = OCr(M) 

15  1F(L)  ly,  19,  16 

16  00  18  1=1, L 

17  IM  = I + M 

D 18  0 = 0 + DCTO( IM,  I ) 

D 19  RC2  = Q 
0 20  RC3R  = UCKK) 

0 21  RNICR  = BPA**K1 

0 22  RN2CR  = S I GM* ( 8NA» *K 1 ) 

D 23  RNINR  = FK2*FK1 *BPASQ 
D 24  RN2NR  = FK2*FK 1*BNASQ 
D 25  RNIR  = RNINR-FLLB 
0 26  RN2R  = RN2NK-FLL8 

C FUD  IS  THE  PRODUCT  ONLY  FOR  RUNNING  INDEX,  HERE  START 

FROM  RC3R. 

D 27  FUD  = RC3R*{  RN1CR*RN1R  + RN2CR*RN2R  ) 

28  IF(M)  39,  39,  29 

29  DO  38  1=1, M 
D 30  FI  = I 

D 31  RC3R  = RC3R*B*(FI+FK)/FI 
D 32  RNICR  = RN1CR*BPA 
D 33  RN2CR  = - RN2CR«BNA 
D 34  RNINR  = { FK2 +F I ) * I FK 1 + F I ) *BPASQ 
D 35  RN2NR  = I FK2 +F I ) * ( FKl+F I ) *BNASQ 
D 36  RNIR  = RNINR  - FLLB 

D 37  RN2R  = RN2NR  - FLLB 

C FUD  IS  running  from  RC3R 

D 38  FUD  = FUD  + RC3R*(  RN1CR*RN1R  + RN2CR*RN2R  ) 

D 39  R = RC1*RC2»FUD 
40  IF(L)  41,  41,  43 
D 41  S = 0. 

42  GO  TO  65 

D 43  SCI  = DCT(L)/(B*»LM2) 

D 44  SC2R  = DCTO(KM,M)*{B**M) 

D 45  SNICR  = BPA**KMi 
D 46  SN2CR  = BNA»»KM1 
D 47  FUD  = 0. 

48  DO  63  I=M1,LM 
D 49  FI  = I 

D 50  SC2R  = SC2R»B*(FK+FI 1/FI 

D 51  Q = 0. 

52  INM  = I-M 

53  DO  55  J=INM,L 

54  JM  = J + M 

0 55  Q = 0 DCTD(JM.J) 

D 56  SC3R  = Q 

D 57  SNICR  = SN1CR*BPA 

D 58  SN2CR  = - SN2CR*BNA 

D 59  SNINK  = (FK2+FI)*(FK1+FI )*BPASQ 

D 60  SN2NR  = ( FK2+F  I ) * ( FKl+FI  ) *B.NASQ 

D 61  SNIR  = SNINR  - FLLB 

D 62  SN2R  = SN2NR  - FLLB 
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C FUD  IS  RUNNING  FROM  SC2 

D 63  FUD  = FUD  + SC2K*SC3R*(  SN1CR»SN1R 
D 64  S = SC1*FUD 

D 65  TCI  = DCT( LM1)*FL2M4* IBI**LM4) 

0 66  TC2R  = DCT(K) 

D 67  rNlR=  BRA**K1 
D 68  FN2R  = S IGM* ( BNA**K  1 ) 

C FUD  IS  RUNNING  FROM  TC2R 

0 69  FUD  = TC2R»(  TNIR  + TN2R  ) 

70  DO  75  1=1, LMl 
D 71  FI  = I 

D 72  TC2R  = TC2R*B»{FK+FI )/FI 

D 73  TNIR  = fNlR*BPA 

0 74  TN2R  = - TN2R*BnA 

C FUD  IN  RUNNING  FROM  TC2R. 

D 75  FUD  = FUD  + TC2R*(  TNIR  + TN2R  ) 

D 76  T = TC1»FUD 

0 77  U = DCT(KLM2)*BISQ*(  (BPA**KLM3)  + 

78  IF(  M-(M/2)*2  ) 79,  79,  89 
D 79  VCl  = DCT(K)*(AI**K1)*(BI»*LM2)*2. 
D 80  Q = DCT(M) 

81  IF(L)  85,  85,  82 

82  DO  84  1=1, L 

83  IM  = I + M 

D 84  Q = Q + DCTDI IM, I ) 

D 85  VNl  = DCT(L)»Q*(  FLL6-FKKA) 

D 86  VN2  = DCTILMI )*FL2M4*BISQ 
D 87  V = VC1*(  VNI+VN2  ) 

88  GO  rO  90 
D 89  V = 0. 

0 90  DHH  = R + S-  T-  U+  V 

D 91  RETURN 
END 


+ SN2CR*SN2R  ) 


SIGL*(BNA**KLM3J  ) 
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FUNCTION  DCT(M) 

C FUNLriGf'l  MAKING  FACfURIALS. 

A IF(M-l)  5,  5,  7 
0 5 OCT  = 1. 

D 6 RLTURN 

7 K = I 

8 IF(M-8)  9,  9,  13 

9 OQ  10  1=2, M 
10  K = K*l 

D ll  DCr  = K 
0 12  RETURN 

13  OU  lA  1=2,8 

14  K=K*I 

0 15  S = K 

16  00  IB  1=9, M 
0 17  FI  = I 

D 18  S = S*FI 
D 19  OCT  = S 
0 20  RETURN 

ENO 


FUNCTION  OCTO(M,N) 

C FUNCTION  MAKING  M FACTORIAL  DEVIDED  BY  N 

C M MUST  BE  LARGER  OR  EQUAL  TO  N. 

C M AND  N MUST  BE  PUAITIVE  INCLUDING  ZERO. 

10  N1  = N+i 

11  IF(M-N)  14,  12,  14 

0 12  DCTD  = 1. 

0 13  RETURN 

14  IF(M)  15,  12,  15 
D 15  S = 1. 

16  DO  18  I=N1,  M 
0 17  FI  = I 

0 18  S = S*FI 

0 19  DCTD  = S 

0 20  RETURN 

ENO 


FACTORIAL. 


APPENDIX  II 


Lower  bounds  using  the  Bazley  Space 


It  was  seen  in  Chapter  V that  for  V > 0,  one  obtains 
t"^  = V“^  - T_ 


(A2.1) 


If  an  upper  bound  satisfies  the  relation 

o 


E°  < £ < E° 
p o p+1  ’ 


(A2.2) 


th 


where  E^  is  the  p eigenvalue  of  from  below,  then  becomes 

positive  definite  if  one  uses  the  pro jection  operator  defined  in  (5.76); 

g 


(A2.3) 


where  g 2p  and  is  an  orthonormalized  eigenfunction  of  with 

the  eigenvalue  E?.  If  the  Bazley  space  is  chosen  by 


Ik  = ...  , 

then  the  corresponding  form  of  (7«ll)  reads 


(A2.h) 


or 


A. . 

ij 


A.  . 

3-J 


< > 
1 J 


1 j 


_lil. 


if  i > g 


if  i ^ g- 


(A2.5) 


(A2.6) 


This  choice  of  the  inner  projection  has  an  advantage  in  that  it  simpli- 
fies the  calculations.  For  this  case 


.-1 


) 
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(A2.7) 
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ki  ^k^  • 


(A2.8) 


The  lowest  eigenvalue  of  the  matrix  obtained  from  (A2.8)  gives  a lower 
bound  to  the  ground  state  energy  of  . 

As  a particular  case,  if  we  let  n ^g,  then  every  element  of  the 
A matrix  has  the  form  (A2.6)  and  the  corresponding  lower  bound  is 
independent  of  an  upper  bound 

The  simplest  example  of  this  choice  is  rendered  by  the  one- 
dimensional Bazley  space 


= t: 


(A2.9) 


then 


(A2.10) 


^^ii^  f^i^  cTkicT^i  ' 


(A2.11) 


A lower  bound,  for  this  case,  is  given  by 


,0/„-l|,„0  V -1  „o 

or  Eg 


E°  + 


which  ever  is  lower. 

For  the  two  electron  ions,  one  has 


i:°  + < lv'^1  T°^“^  = -Z^+ 
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and 


2 5 2 
" ■ 8^ 


It  is  seen  that  for  Z = 1,  E°  is  a lower  bound  to  the  ground  state 

2 l6 

energy  whereas  for  Z 212,  - Z + ^ lower  bound.  It  is 

shown*-^^  that  + < I'®!  v"^l  1'°  is  a lower  bound  to  the 

ground  state  energy  of  if 


^l<  So  <4 
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